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Abstract: We use gauge-gravity duality to compute spectral functions of fermionic operators 
in a strongly-coupled defect field theory in p-wave superfluid states. The field theory is (3-1-1)- 
dimensional = 4 supersymmetric SU{Nc) Yang-Mills theory, in the 't Hooft limit and with 
large coupling, coupled to two massless flavors of (2-|-l)-dimensional A/" = 4 supersymmetric 
matter. We show that a sufficiently large chemical potential for a U{1) subgroup of the 
global SU (2) isospin symmetry triggers a phase transition to a p-wave superfluid state, and 
in that state we compute spectral functions for the fermionic superpartners of mesons valued 
in the adjoint of SU{2) isospin. In the spectral functions we see the breaking of rotational 
symmetry and the emergence of a Fermi surface comprised of isolated points as we cool the 
system through the superfluid phase transition. The dual gravitational description is two 
coincident probe D5-branes in AdS^ x with non-trivial worldvolume SU{2) gauge fields. 
We extract spectral functions from solutions of the linearized equations of motion for the D5- 
branes' worldvolume fermions, which couple to one another through the worldvolume gauge 
field. We develop an efficient method to compute retarded Green's functions from a system of 
coupled bulk fermions. We also perform the holographic renormalization of free bulk fermions 
in any asymptotically Euclidean AdS space. 
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1. Introduction and Summary 

The Anti-de Sitter /Conformal Field Theory correspondence (AdS/CFT) ||^, |2|, ^, and more 
generally gauge-gravity duality, is a holographic duality between a weakly-coupled theory of 
gravity in some spacetime and a strongly-coupled field theory living on the boundary of that 
spacetime. Gauge-gravity duality thus provides a powerful new tool for studying strongly- 
coupled, scale-invariant field theories in states with finite charge density, and hence may be 
useful in condensed matter physics, for instance in understanding low-temperature systems 
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near quantum criticality |5|, ^, ^. In particular, many special properties of certain high- 
Tc superconducting materials may be due to an underlying quantum critical point [§, |5|, ^. 
Gauge-gravity duality may provide valuable insight into the physics of such materials. 

Of central importance for potential condensed matter applications is the holographic 
description of a Fermi surface^ P, 10, lT| T3]. On the field theory side, the minimal ingredients 



are some strongly-coupled theory with a global U{1) symmetry, in a zero-temperature state 
with a finite U{1) chemical potential, and some fermionic operator charged under the U{1). 
Holographic calculations of the fermionic spectral function, as a function of frequency and 
momentum, reveal a pole at zero frequency but finite momentum, which defines the Fermi 
momentum. The pole represents an excitation about a Fermi surface. 

On the gravity side, the minimal ingredients are gravity and a U{1) gauge field, plus some 
bulk Dirac fermion charged under the U{1). The bulk geometry is a Reissner- Nordstrom black 
hole. The bulk fermion is dual to the fermionic operator, and the spectral function of the 
operator is extracted from solutions of the linearized bulk equation of motion, the Dirac 
equation. These Fermi liquids are, generically, not Landau Fermi liquids, although the exact 
properties depend on the mass and charge of the bulk fermion. 

The bulk theory can also describe a phase transition to s-wave superfluid states, if a scalar 
charged under the U{1) is present |13, On the gravity side, the Reissner-Nordstrom 



black hole grows scalar hair at low temperature, that is, a solution with a non-trivial scalar 
becomes thermodynamically preferred to Reissner-Nordstrom. In the dual field theory, the 
thermodynamically-preferred state includes a nonzero expectation value for a scalar operator 
charged under the global U{1), which we will refer to as the operator "condensing." The 



phase transition is second order with mean-field exponents |14, 15|. 

Gauge-gravity duality can also describe p-wave superfluids, that is, superfluids in which 
the condensing operator is a vector charged under the U{1), thus breaking not only the U{1) 
but also rotational symmetry (to some subgroup) On the gravity side, the minimal 

ingredients are gravity and non-Abelian gauge fields. The simplest case is an AdS geometry 
and SU (2) gauge fields, A^j, with Lorentz index M and a = 1, 2, 3 labels the SU{2) generators 
Ta- Here the U{1) is a subgroup of SU{2), for example the U{1) in the direction, which we 
will call 17(1)3. At high temperature the thermodynamically preferred geometry is Reissner- 
Nordstrom with nonzero A^. At low temperature, the charged black hole grows vector hair: 
the preferred solution has non-trivial A],. The dual field theory has three conserved currents, 
Ja, dual to the gauge fields. A chemical potential, producing a finite density (J3), explicitly 
breaks SU{2) to C/(l)3, and the transition occurs at large chemical potential, where the 
thermodynamically preferred state has nonzero (Jf ). 

In bulk calculations for both the s- and p-wave, a major technical simplification is the so- 
called probe limit, in which the charge of the bulk scalar, or the SU{2) Yang-Mills coupling, is 
sent to infinity, so that the scalar or Yang-Mills stress-energy tensor on the right-hand side of 
Einstein's equation becomes negligible. The bulk calculation then reduces to solving the scalar 



^For an alternative approach, see ref. 
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or Yang-Mills equation of motion in a fixed Reissner-Nordstrom background. The probe limit 
is sufficient to detect the transitions and determine that they are second order. In either case, 
however, if we cool the system then, as shown in refs. [^], the matter fields' stress-energy 
tensor grows and we can no longer trust the probe limit. Reaching zero temperature requires 
solving the fully-coupled equations, as done in refs. |]l5|, 17, |l8|. 



The zero-temperature limits of the bulk hairy black hole solutions generically involve 
a domain wall interpolating between two regions, one near the boundary and one deep in 
the interior of the spacetime. For example, the geometry may interpolate between a near- 
boundary AdS space and an interior AdS space with a different radius of curvature and speed 
of light [19, 2C, 22, 23]. In field theory language, the interior AdS space represents an 
emergent conformal symmetry at low temperature and finite charge density. In other words, 
the emergent AdS represents a quantum critical point. 

Holographic calculations of fermionic spectral functions in zero-temperature s-wave su- 
perfluid states |24, |2^, 26, 27 1 exhibit the so-called 'peak-dip- hump' structure pi], expected 



to be relevant in high- Tc superconductors [28|, as well as, for suitable mass and charge of the 
bulk fermion, continuous bands of poles |26] and, for suitable coupling to the bulk scalar, a 



gap, i.e. poles in the spectral function at nonzero momentum and nonzero frequency |2c 

Generally, the bulk actions used in holographic constructions of superfluids and Fermi 
surfaces are not derived from any particular string theory construction. In other words, 
they are basically ad hoc models built from the minimal ingredients needed to capture the 
essential physics. Simple models have one big advantage (besides simplicity!), namely a kind 
of universality: the results may be the same for many different theories, regardless of the 
details of their dynamics. 

On the other hand, knowing the detailed dynamics of a specific dual theory, meaning 
the fundamental fields and Lagrangian of some microscopic, weak-coupling description, also 
has advantages. For example, the holographic results may tell us that a superfluid phase 
transition occurs, but may not tell us why. Is a nonzero (Jf ) the result of some pairing 
mechanism? If so, is the pairing mechanism the same in every dual theory? Knowing an 
exact dual theory may help to answer such questions, for example by providing some weak- 
coupling intuition.-^ Finding a dual Lagrangian means "embedding" the bulk theory into a full 
string or supergravity construction, built for example from D-branes (for which we know the 
worldvolume theories). String (and M-) theory embeddings of holographic s-wave superfluids 
appear in refs. ||, |o|, |3l|. 

A string theory embedding of holographic p-wave superfluids, in the probe limit, appears 
in refs. |32, |3^, 34, |3^]. Here we begin with Nc Dq-branes and Nf Dp-branes. Taking 
the usual decoupling limit for the Dq-branes, which in particular means Nc — )• oo, we obtain 
supergravity in the near-horizon geometry of the Dq-branes. Non-extremal Dq-branes produce 
a black hole geometry. The probe limit consists of keeping Nf flxed as Nc — )• oo, such 
that Nf/Nc — )• 0. The dynamics of Nf coincident Dp-branes is then described by the non- 



good recent example is ref. |2E|]. 
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Abelian Born-Infeld action (plus Wess-Zumino terms) in the near-horizon Dq-brane geometry. 
Truncating that action to leading order in the field strength, we obtain a Yang-Mills action 
in a black hole geometry. 

We then know precisely what the dual field theory is: the Dq-brane worldvolume theory, 
with gauge group SU{Nc), in the large- A^c and strong coupling limits, coupled to a number 
Nf of fields in the fundamental representation of the gauge group, i.e. flavor fields. We 
will call this the Dq/Dp theory. If the Dp-branes do not overlap with all q spatial Dq-brane 
directions, then the flavor fields will only propagate along a defect. The probe limit consists 
of neglecting quantum effects due to the flavor fields, such as the running of the coupling, 
because these are suppressed by powers of Nf/N^. These theories generically have bound 
states similar to mesons in Quantum Chromodynamics (QCD). The U{Nf) gauge invariance 
of the Dp-branes is dual to a global U{Nf) analogous to the isospin symmetry of QCD. In 
such systems the p-wave transition occurs when a sufficiently large isospin chemical potential 
triggers vector meson condensation (as we review in section ^). 

Our goal is to use such a string theory system to compute fermionic spectral functions 
in the p-wave phase. 

We choose our Dq-branes to be D3-branes. The dual field theory is then a CFT, A/" = 4 
supersymmetric Yang- Mills (SYM) theory with gauge group SU (Ac) in the 't Hooft limit of 
Ac — )• oo with large 't Hooft coupling, A = gY]\,jNc — )• oo. The near-horizon gravity solution 
is (4-|-l)-dimensional AdS times a five-sphere, AdS^ x S^, with Ac units of Ramond-Ramond 
(RR) five- form flux on the S^. At finite temperature AdS becomes AdS-Schwarzschild. 

We consider supersymetric probe Dp-branes extended along AdSp x , where super- 
symmetry requires |P — Q| = 2 |36|. We focus on P > 3, since only in those cases is a vector 
condensate ( Jf ) possible.^ We will study only trivial embeddings of such Dp-branes, that is, 
we will only study solutions in which all Dp-brane worldvolume scalars are trivial. 

Our bulk fermions will be the Dp-branes' worldvolume fermions. These fermions are in 
a supermultiplet with the worldvolume scalars and gauge field, hence they are in the adjoint 
of the worldvolume U{Nf), and couple to the gauge field via the gauge-covariant derivative.^ 
In other words, supersymmetry determines the charges of the fermions. For example, we will 
use A^^ = 2, where we find three fermions with charges +1, —1 and under [7(1)3. 

To compute fermionic spectral functions we need the linearized equations of motion, the 
Dirac equation, for these fermions. Fortunately, the fermionic part of the D-brane action, 
for D-branes in arbitrary backgrounds (including RR fields) is known to quadratic order 
|39t| . The form of the action is determined by supersymmetry and T-duality |3^] , as 
we review in section |. For our Dp-branes extended along AdSp x S^, we perform a reduction 
on the to obtain a Dirac equation in AdSp, following ref. very closely. The spectrum 
of AdSp fermion masses are fixed by P, Q and the coupling to the background RR five-form. 

We emphasize a major difference between our systems and the models of refs. IC, [Tl[ |: 
in our embedding of the Dirac equation into string theory, the mass and charge of the fermions 

^One exception is a D5-brane along AdS2 x S'^, which we study in section 4.1 (but not in p-wave states). 
''Like all worldvolume fields, they are not charged under the diagonal (7(1) C U{Nf). 
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are fixed by super symmetry and T-duality. We are not free to dial the values of the mass and 
charge, unlike refs. |lO|, pl|] . 

Much of our analysis will be valid for any supersymmetric Dp-brane extended along 
AdSp X S^, with P > 3, but one particular Dp-brane is attractive for a number of reasons, 
namely the D5-brane extended along AdS^ x S'^ (P = 4 and Q = 2). From the bulk point of 
view, this D5-brane is the only Dp-brane with a massless worldvolume fermion, as we show 
in section 4.1. That makes both our numerical analysis, and comparison to refs. |ll[ (in 
which the fermions were massless), much easier. 

With two coincident P = 4, Q = 2 D5-branes, the dual field theory is (3+l)-dimensional 
= 4 SYM coupled to Nf = 2 massless (2-|-l)-dimensional = 4 supersymmetric flavor 
fields. The classical Lagrangian of the theory, with couplings that preserve the SO{3,2) 
conformal symmetry of the (2-|-l)-dimensional defect, appears explicitly in refs. |4^ . 
We write the explicit form of the fermionic operators dual to the D5-branes' worldvolume 



fermions in section |4.2| , following refs. [40, 41 1 very closely. These fermionic operators are 
mesinos, the supersymmetric partners of mesons. 

The P = 4, Q = 2 D5-brane is also attractive for potential condensed matter appli- 
cations. As mentioned in ref. |43|, many real condensed matter systems are effectively 
(2-|-l)-dimensional degrees of freedom interacting with ambient (3-1-1 )-dimensional degrees 
of freedom. The D3/D5 theory also exhibits a rich phase structure, explored in detail in 
refs. m, H, m, m, H^, H, H, l^, including for example a Berezinskii-Kosterlitz-Thouless 
transition (with finite charge density and magnetic field for the diagonal U{1) C U{2)) |5l[| . 
Here we focus on the D3/D5 theory's phase diagram with finite isospin chemical potential. 

As always in the probe limit, we cannot access the T = 0, finite chemical potential 
ground state. The P = 4, Q = 2 D5-brane is again attractive, however, because we know that, 
unlike many Dp-branes, with zero temperature and zero chemical potential, fully back-reacted 



solutions appear to preserve an AdS factor in the geometry, namely an AdS4 p2, 53]. That 



suggests that the field theory retains 50(3,2) conformal invariance even including quantum 
effects due to the fiavor, which was indeed proven in ref. |42|. Whether some scale invariance 
emerges with zero temperature and finite chemical potential is unclear. 

On a technical level, our goal is to solve the Dirac equation for a massless fermion in the 
adjoint of SU{2) confined to an AdSp submanifold of (4-|-l)-dimensional AdS-Schwarzschild. 
For any Dp-brane, the three worldvolume fermions decouple in the normal (non-superfluid) 
phase, where Al. is zero, but couple to one another in the superfluid phase, where A], is 
nonzero. These couplings indicate that, in the field theory, the dual fermionic operators 



experience operator mixing under renormalization group flow 55 1. In the fleld theory, 
the retarded Green's function, and hence the spectral function, becomes a matrix with off- 
diagonal entries. 

We thus develop a method to compute the retarded Green's function for bulk fermions 
coupled to one another. Our method is essentially a combination of the method of ref. [54, 55[, 
for coupled bosonic flelds, with the method of ref. [^, for free fermions. Our method 
is actually very general, i.e. applicable to any system of coupled bulk fermions, not just to 
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fermions on the worldvolume of probe Dp-branes, and is especially convenient for numerical 
analysis. We thus explain our method first, in section]^. 

As an added bonus, we also perform, to our knowledge for the first time, holographic 
renormalization for fermions in AdS.^ More precisely, we study a single free fermion in any 
space that asymptotically approaches Euclidean-signature AdS and determine the countert- 
erms needed to render the on-shell action finite without spoiling the stationarity of the action. 
Our results rigorously justify many of the ad hoc prescriptions used in the literature, where 
divergences of the on-shell action were simply discarded. 

For the P = 4, Q = 2 D5-brane, using our method for coupled bulk fermions, we 
numerically compute spectral functions for mesinos as we cool the system through the p- 
wave superfluid phase transition. Due to the operator mixing, or equivalently the coupling of 
the fermions in the bulk, we see that the spectral function of even a neutral fermion develops 
a nontrivial feature, a peak, as the system enters the p-wave phase. 

Furthermore, as we lower the temperature, the zero- frequency spectral measure^ is clearly 
no longer rotationally invariant, and in fact at the lowest temperatures we can reliably access 
in the probe limit, the main features of the spectral measure are five largely isolated peaks 
in the {kx,ky) plane, two on the kx axis, two on the ky axis, and one at the origin. These 



results are very similar to the T = results of ref. |18], where the bulk theory was gravity 
and SU{2) gauge fields in (3-|-l)-dimensions, in the T = vector-hairy black hole geometry. 
In that case, for a fermion in the fundamental representation of SU{2), the spectral measure 
consisted of two points on the kx axis, located symmetrically about the origin. The prediction 



of ref. 1 15] for fermions in the adjoint representation would be three points on the kx axis, 
one at the origin and two at finite kx, positioned symmetrically about the origin. At finite 
temperature we see five points emerging, but we strongly suspect that, if we could access the 
T = limit, we would indeed see only three points, as we discuss in section ^ 

We cannot resist drawing an analogy between our system and certain experimentally- 
realized p-wave superconductors (see also ref. [^]).'' In that context, a "reduction of the 
Fermi surface" to certain points in momentum space has been proposed for the ruthenate 



compound Sr2Ru04, |59|: the p-wave state is supported by ferromagnetic fluctuations that 
increase the propensity for electrons to form spin triplet Cooper pairs, with an odd (p-wave) 
Cooper pair wave function.^ Scattering channels with momentum transfer Q = (0, 0), as is the 
case in a ferromagnet, should be enhanced in the system, as opposed to scattering channels 
oi Q = (vTjTr), which is the case in an anti- ferromagnet. Small momentum transfer is best 
accomplished by a strongly peaked density of states at the Fermi level, as occurs for example 
with van Hove singularities, where the density of states diverges. This lies at the heart of the 



^For the holographic renormahzation of fermions in Schrodinger spacetime, see ref. | |57| . 

® As mentioned above, the retarded Green's function, and hence the spectral function, is genericaUy a matrix. 
The spectral measure is simply the trace of the spectral function. 

'^We thank Ronny Thomale for many useful conversations about real p-wave superconductors. 

*This is rather particular, bearing in mind that a large number of generic spin interactions, for example 
induced by superexchange processes, favor antiferromagnetic fluctuations. 
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strong suspicion that a Fermi surface localized to certain points with a high density of states 
may account for a suitable setup to support p-wave pairing. 

The paper is organised as follows. In section ^, we describe our method for computing 
retarded Green's functions for coupled bulk fermions. In section ^ we review general features 
of Dq/Dp holographic p-wave superfluids and demonstrate a p-wave transition using the 
P = 4, Q = 2 D5-brane. In section ^ we write the fermionic part of the Dp-brane action, 
perform the reduction of the worldvolume Dirac equation to AdSp, and, for the P = 4, 
Q = 2 D5-brane, match bulk fermions to dual field theory operators. In section ^ we present 
our numerical results for the fermionic retarded Green's functions using the P = 4, Q = 2 
D5-brane. We conclude with suggestions for future research in section |^. The holographic 
renormalization of fermions in AdS appears in the appendix. 

Section |2| and the appendix are technical, and may be read independently of the rest 
of the paper. Readers who only want to understand our system and our numerical results 
should read at least sections 3.2, 4.3 and S. 



2. Holographic Fermionic Operator Mixing 
2.1 Revievi^: Free Fermions 

We begin by studying a single free fermion in AdS space. In particular we will review how to 
extract the field theory fermionic two-point function from a solution for a bulk Dirac fermion, 
following refs. |§. 

In this section we will work mainly with Euclidean-signature AdS space, with the metric 
written in Fefferman-Graham form,^ 

du^ 1 

ds"^ = 9AB dx^dx^ = — TT H — o^a dx^dxK (2.1) 

The boundary is at u = 0. Notice that throughout the paper we use units in which the radius 
of AdSd+i is equal to one. 

We will study a bulk Dirac spinor ^ . The Dirac action (plus boundary terms) is 

S = j d'^+^x ^{^p^-m ^>^>) + Sbdy, (2.2) 

where, picking one of the spatial directions to be "time," with corresponding 7*, we define 
^ = ^^j^. We write the AdSd+i Dirac operator P below. Here Sbdy includes boundary terms 
that do not affect the equation of motion. 

The AdS/CFT correspondence is the statement that a theory of dynamical gravity on 
AdSd+i is equivalent to a d (spacetime) dimensional CFT that "lives" on the boundary of 
AdSd+i- Every bulk field is dual to some operator in the boundary CFT. The precise state- 
ment of the correspondence equates the on-shell bulk action with the generating functional 

^Capital Latin letters A,B,... will always denote all the AdSd+i directions, including the radial direction 
u, while lower-case Latin letters will denote field theory directions: i, j = 1, . . . d. 
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of connected CFT correlation functions. The bulk field \I' is dual to some fermionic operator 
O in the dual d-dimensional field theory. The on-shell bulk action, S, acts as the generating 
functional for correlators involving O. In other words, to compute renormalized correlators 
of O, we take functional derivatives of S with respect to some source. 

Generically, however, both the on-shell bulk action and the CFT generating functional 
diverge. On the bulk side, the divergences arise from the infinite volume of AdSd+i, i-e. they 
are long-distance or infrared (IR) divergences. In the field theory the divergences are short- 
distance, ultraviolet (UV) divergences. To make the AdS/CFT correspondence meaningful 
we must regulate and renormalize these divergences. 

Holographic renormalization proceeds as follows (see ref. |61] and references therein). We 



first regulate the on-shell bulk action by introducing a cutoff on the integration in the radial 
direction: we integrate not to it = but to some u = e. We then add counterterms on the 
u = e surface to cancel any terms that diverge as we remove the regulator by taking e — )• 0. 
Generically, the form of the counterterms is fixed by symmetries, and the coefficients of the 
counterterms are adjusted to cancel the divergences. Once the counterterms are known, we 
can proceed to compute functional derivatives of the on-shell bulk action, always taking e — )• 
in the end, thus obtaining renormalized CFT correlation functions in a way that is manifestly 
covariant and preserves all symmetries. 



As first observed in ref. [62|, when we evaluate the Dirac action on a solution, the bulk 
term obviously vanishes. The nonzero contribution to the on-shell action comes from Sbdy, 
which involves terms localized on the u = e surface. As observed in ref. ||6^, the form of S^y 
is fixed by demanding a well-defined variational principle for the Dirac action. Formally, S^ar 
thus includes two types of terms, 

Shdy = Syar + ScT, (2.3) 



where S^ar are the terms required for the variation of the action to be well-defined |63], while 
ScT are the counterterms, which do not affect the variation of the action. 

In the appendix we perform the holographic renormalization of the Dirac action. In 
particular, we determine the counterterms in Sct- The details of holographic renormaliza- 



tion are well-known for various species of bulk fields, for example for the metric [34|, scalar 
fields [0], and gauge fields ||6^. To our knowledge, the only detailed analysis of holographic 
renormalization for fermions was in the (more complicated) context of non-relativistic gauge- 



gravity duality, in ref. |57]. As shown in the appendix, however, in the relativistic case the 
holographic renormalization procedure for fermions very closely parallels the procedure for 
scalar s. 

As shown in the appendix, the details of the holographic renormalization depend on the 
value of m. Some values of m are special, for example when m is half-integer (in units of 
the AdS(i+i radius), counterterms logarithmic in e (rather than just polynomial in e) are 
needed. For simplicity, in this section we will restrict to values of m that are positive and not 
half-integer. Our arguments are easy to generalize to any value of m. 
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In this section we will also restrict to four- and five-dimensional AdS spaces, which we 
will collectively denote as AdS^+i with d = 3,4, primarily for pedagogical reasons: in these 
cases the bulk Dirac spinor has four complex components, and we can write explicit 4x4 bulk 
Dirac F-matrices. Additionally, we note that AdSd+i spaces with d < 4 are the cases most 
relevant for condensed matter applications (as opposed to, say, AdSj). The generalization to 
other dimensions is straightforward. In the appendix we work with arbitrary d. 

In later sections we will be interested in computing finite-temperature, real-time correla- 
tion functions, in particular the retarded Green's functions, in which case the bulk geometry 
will be Lorentzian-signature AdS-Schwarzschild. We review the prescription for obtaining the 
retarded Green's function in such cases at the end of this subsection. 

Varying the above action we obtain the bulk equation of motion, the Dirac equation. 



M 



^ - = 0, 



(2.4) 



where = u6^ are the inverse vielbeins associated with the metric in eq. (2.1).^^^ The 



curved-space covariant derivative is 

Dm = 



(2.5) 



where {u}m)ab ^P^'^ connection associated with the metric in eq. (p.l|). The only 

nonzero components of the spin connection are (wj)^^ = ^ 6ij, so that Du = du and the other 
components of Dm are 

1 T 

(2.6) 



D, = di + l- [7",y 



4m 



We can now simplify the Dirac equation 





e^j^DM^ 



m ^ 



M 



U7 



n7 



om 



4 

2 ' 



m 



(2.7) 



We will work with a single Fourier mode, so we let ^' — )• e ' where, without loss of 
generality, we have chosen the momentum to point in the x direction. The Dirac equation 
is then 

ui^du + iku 7^' 7" 



m 



^ = 0. 



(2. 



We will now choose an explicit basis for the F-matrices. We will use a basis in which all 
the F-matrices are Hermitian, 



7 




7 




7 




(2.9) 



^"Recall that for inverse vielbeins, the upper index is general coordinate and the lower index is local Lorentz. 
The 7^ obey the usual algebra {7'*,7^} = 2 5^^^. 

In a p-wave superfiuid phase rotational symmetry is broken, so there, to study the most general case, we 
must use a momentum with nonzero components in different directions, as we will discuss in section ^ 
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where ai, o"2 and as are the usual Pauh matrices, 

f l\ f -i 



cri 



0-2 



10/ \i I 

Next we wih define two sets of projectors. The first set is 

/O \ 



1 
-1 



(2.10) 



n+ = -(i + 7") 



1 







n- = ^(i-7") 



V V 



/I \ 



1 

V 0/ 



(2.11) 



We use these to define = ^ (1 it 7") ^ so that 7"^'^ = it^'-i-. The second set of projectors 
was used for example in refs. 0, O], 



/O 



Hi 



1 



(1 + ^7-7*7-) 



\ 







1 



/I 



(i-^Vt") 



1 







(2.12) 



To make converting between ^'j- and ^'1^2 easy, we explicitly write ^ first as ^'_|- + 'I'- and 
then as ^1 + ^2, 

/ \ /^-u\ / \ /^2u\ 







+ 





\ J 







+ 





V / 



(2.13) 



where the subscripts u and d indicate the "up" and "down" components of the effectively 
two-component ^± and ^'1,2- Identifications such as "f+u = 'I'2d are then obvious. 

We have a choice of whether to use ^± or ^1,21 although of course, we can easily translate 
between the two options using eq. ( 2.13| ). We will choose whatever is most convenient for a 
given question. 

For example, the projectors ni^2 commute with the operator in eq. (|2.8| ), which tells us 
that, for a free fermion, the equations for ^1^2 decouple. That makes ^^1^2 especially attractive 
for numerical analysis, hence we employ them in sections ^ and ^.^^ Explicitly, the equations 
for ^'i 2 are 



udu- -+ma3 



ku 



udu 



+ mas + ku 



^-1 



0, 



0. 



(2.14) 
(2.15) 



^■^As mentioned in footnote |l^, in the p-wave superfluid phase, the most general momentum has nonzero 
components in multiple directions. That means ^/i and ^2 will no longer decouple because other F-matrices, 
such as 7'', will appear in the equation of motion, and these do not commute with Hi, 2- Nevertheless, when 
studying the p-wave superfluid phase we use »I'i,2 to make the comparison with the rotationally-symmetric 
case easier. 
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On the other hand, the asymptotic behavior of ^ is most succinctly described using 
hence we use these frequently below, especially in the appendix. In terms of "i>±, the equation 
of motion becomes 

udu-^-m^'^+ + kua3'^^ = 0, (2.16) 
udu-^ + m^'^- + kua3'^+ = 0. (2.17) 
These first-order equations give rise to the second-order equations 



u \ 4 2' 

The leading asymptotic behaviors of ^± are 

*± = c±(A;)ni±™ + 0Li+i±") . (2.19) 



= 0. (2.18) 



where c±{k) are spinors that obey Ii±c±{k) = ±c±{k), and which may depend on k, as 
indicated. 

As reviewed above, to compute renormalized correlators of the dual operator O, we take 
functional derivatives of S with respect to some source. We identify the source for O as the 
coefficient of the dominant term in ^''s near-boundary expansion (the term that grows most 
quickly as n — )• 0). From eq. ( ^) , we see that the dominant term is the U2 ™ term, hence 
we identify c-{k) as the source for O. More formally, we equate 



g-5ren[c_,C_] ^ ^^^p 



(2.20) 



where the left-hand-side is the exponential of minus the action in eq. ( p. 21) , evaluated on 
a solution and properly renormalized (hence the subscript), and the right-hand-side is the 
generating functional of the dual field theory, with c_(/c) acting as the source for the operator 
0.13 Upon taking minus the logarithm of both sides, we find that the on-shell bulk action is 
the generator of connected correlators. 

For bulk bosonic fields, we must solve a straightforward Dirichlet problem: we fix the 
leading asymptotic value of the field, allow the field to vary, and then impose a regularity 
condition in the interior of the space to fix the entire solution. This procedure is dual to 
the statement that once we choose a source, the dynamics of the theory determines the 
expectation values of the dual operator. 

The story for fermions is more subtle, because ^'+(n, k) and ^-(n, k) are not independent 



63 1 . Each one determines the canonical momentum associated with the other (see for 



example ref. ]6^). In the bulk Dirichlet problem, then, we cannot fix their asymptotic values 



^^As we review in the appendix, for a bulk fermion with mass m, in a standard quantization the dimension 
AofOisA = | + In the appendix we also discuss the chirality of O (when d is even). 
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c±{k) simultaneously, but can fix only one, the coefficient of the dominant term, C-{k), and 
then vary the field. As shown in refs. [^], for the action to remain stationary under such 
variations, we must add a boundary term to the action. 



Svar = J d'^X^^ + ^^, (2.21) 



where the integration is over the u = e hypersurface, ^77 = e~'^ is the square root of the 
determinant of the induced metric at u = e, and ^± are evaluated at u = e. 

Indeed, since the bulk action is ffi'st-order in derivatives, the only nonzero contribution 
to the on-shell action comes from the boundary terms Shdy = S^jar + Sct- Generically, 
when evaluated on a solution, divergent terms appear in Syar, which are canceled by the 
counterterms in Sct- Notice that, to preserve stationarity of the action, Sct must involve 
only ^'-(e, k), since that is held fixed under variations. We write the counterterms explicitly 
in the appendix. 

The principal result of the appendix is the renormalized on-shell action: we evaluate Sbdy 
on a solution and take e — t- to obtain (for positive, non-half-integer m) 

-id 



d'^xc+c-, (2.22) 

We can now easily compute the renormalized connected correlators of O and O by taking 
functional derivatives of Sren- For example, the renormalized one-point function of O is 

If we use the fact that the on-shell bulk action must be Hermitian, S = S\ then we also have 

Sren = Sl^n = j ^'^^ [c+C-]^ = j d'^XC^C+, 

hence we also find, as we should, 

{0)ren = = "C^- (2-24) 

We can obtain two-point functions via second functional derivatives, for example 

The equation of motion plus some regularity condition in the interior of the spacetime will 
relate c+ and c_ (recalling that we fix c_ and vary c+). The equation is linear, hence the 
relation will be linear: c+ = —G{k) 7* c_, for some matrix G{k) which will turn out to be the 
Euclidean Green's function. We include a factor of 7* because, as discussed in refs. [^,11], 
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the Euclidean Green's function is actually (O C^)^g„) which differs from (^O 
of 7*. We indeed find 



(OO) 



G{k). 



by a factor 



(2.26) 



In general, we must extract G{k) 7* from a solution by imposing some regularity condition in 
the bulk of the spacetime (in our coordinates, the u — ?■ 00 region), which fixes Cj^ in terms 
of c_. We review that procedure for Euclidean AdSd+i in the appendix and for Lorentzian- 
signature AdS-Schwarzschild below. 

We can also reproduce the formulas used in refs. ^ by switching to 'I'i,2- In that 
case, the equations for ^1 and ^'2 decouple, hence in the Green's function the Hi and 112 
subspaces will not mix. Writing cj^ = —G{k) 7* c_ explicitly, we will have (suppressing the k 
dependence of c±{k)) 

/ \ /O 1 \ 





V+dJ 



G22{k) I2 



^11(^)12, 



1 



V 



1 
10/ 





V J 



G22{k)l2 



Gii{k) I2 



/ \ 



\C-dJ 



where blank entries represent zero, I2 is the 2x2 identity matrix, and Gn and G22 represent 
the components of the Green's function in the Hi and 112 subspaces, respectively. Given a 
bulk solution for we obtain the Green's functions simply by reading off the asymptotic 
values of c+(A;) and C-{k) and then constructing 

G22(fc) = -^, Gn(A;) = -^. (2.27) 

Finally, we review the prescription of ref. |6C] to compute the retarded two-point function 
in the finite-temperature, Lorentzian-signature case. Here the geometry is AdS-Schwarzschild, 
with a horizon at some position Uh- To obtain the retarded two-point function, we require 
that, near the horizon, the bulk solution for ^ has the form of wave traveling into the horizon 
(out of the spacetime), i.e. an in-going wave. The asymptotic form for ^ near the boundary 
is the same as in eq. ( ^.191) (for positive, non-half- integer m). Following ref. [|^], in the 
regime of linear response, we have 

c+{uj,k) = -iG^{u;,k)-f* c-{u},k), (2.28) 

where G^{uj, k) is the retarded Green's function. Notice that here we distinguish the frequency 
uo from the momentum k, and 7* is now anti-Hermitian, 



7 




Eq. ( |2.28| ) is essentially just an analytic continuation from the Euclidean case: 7* 
a free fermion, we obtain (see also eq. (A17) of ref. |pT|]) 

^, Gfi(a;,A:) = ^. 

C-—U C—d 



^22 



[u,k) 



(2.29) 



^7*. For 



(2.30) 
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2.2 Coupled Fermions 

We now consider multiple bulk fermions, say N of them, with a = 1, . . . ,N, coupled to 
one another. The fact that the linearized fluctuation of the couple in the bulk is dual to 
the statement that the fermionic operators in the field theory mix with one another under 
renormalization group flow. 

We will work in Lorentzian signature, and finite temperature, so that the bulk geometry is 
AdS-Schwazrschild, with a horizon at some position Uh- We consider fermions with quadratic 
couplings of the form (with implicit summation over a, b) 

S = i j <f+^X ^ {^a P'^a - ^aKb'^b) + Sbdy, (2.31) 

for some matrix Aab that need not be diagonal in either the a, b indices or the spinor indices. 
As a concrete example, in later sections we will introduce a bulk SU{2) gauge field Am and 
a bulk fermion valued in the adjoint of SU (2). The indices a, b are then SU (2) indices, hence 
we will have three bulk fermions (for ri, T2, and T3) with a coupling, coming from the gauge- 
covariant derivative, of the form Cafec^a e^7^ {Am)i,^c, which is obviously not diagonal in 
either SU{2) indices or in spinor indices (because of the 7"^). 

For the following arguments, we do not need to know any details about the equations of 
motion. We will only exploit one important feature. Using the II-i- projectors, we will always 
obtain equations similar to eqs. ( 2.16| ) and ( 2.17] ). We will then always be able to write these 



equations in the form 

Vab±'^b± = Mac±'^cT, (2.32) 

where Vafe± is some differential operator, involving in particular du, and Mac± is a matrix 
representing the couplings among not only the ^'q, which come from Aab, but also the terms 
from P'^a that produce couplings between '^a+ and ^'a-, for example the terms proportional 
to the momentum k in eqs. (2.16) and ( p.l7| ). The key feature is that only the '^a± are on 



the left-hand-side, while only the ^'^^ are on the right-hand-side. 

In practical terms, the total number of complex functions for which we must solve is 
4 X A^, since each has four complex components. In other words, we need to decompose 
the not only into and ^a-, but also into the up and down components, ^a+u, ^a+d; 
^a-u, and ^a-d- When convenient, we may sometimes think of eq. ( ^.32 ) as equations 



describing these 4 x A^ coupled functions, which we may sometimes refer to as "fields." 

Clearly, if we solve for all the insert the solutions into the bulk action, and take 
functional derivatives, we will obtain field theory retarded Green's functions that are matrices, 

{u),k). In principle, we may be able to diagonalize the equations of motion and obtain 
decoupled equations, in which case the Green's function will be diagonal. Given the bulk 
solutions for the ^'a, we then extract the elements of G^^, (cj, A;) using eq. (2.30). In some cases. 



however, diagonalizing the equations of motion may be prohibitively difficult, i.e. practically 
impossible. We can always resort to numerics to find solutions, but we will then be forced to 
compute elements of the un-diagonalized {uj,k). We thus need to know what combinations 
of the asymptotic values Ca+ and Ca- give an arbitrary element G^^{u},k). 
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We will describe a prescription to obtain the matrix G^^{ui,k), assuming we have bulk 



solutions for the The method is a hybrid of the methods in refs. |]54| , 55 1 and refs. |5i 



Refs. p% pq ] described a general method to construct a retarded Green's function for coupled 
bulk scalar and gauge fields, while refs. [10, 56] described general methods for computing 
Green's functions from fermions in the bulk. 

The first observation is that we can construct second-order equations for the bulk fields, 
the ^a±-, that will be similar to eq. ( p. 18 ). We actually don't care about the exact form of 



these equations. We only need to know that such equations exist. We thus have a system of 
2N second-order linear equations, for which we expect 2 x 2N linearly-independent solutions. 
We must therefore fix two boundary conditions for each field to specify a solution for the 



entire system. Following refs. [55|, we fix these boundary conditions near the horizon Uh- For 



example, the ^a+u will have the near-horizon form 

^'a+« = na+u [u - Uh)"^ + .... (2.33) 

where Ua+u and a are constants (independent of u) and . . . represents terms that decay faster, 
as n — )• u/i, than the terms shown. The two constants Ua+u and a are the two degrees of 
freedom we have to specify the solution. Generically, the equation of motion will only be 
satisfied for two values of a, one describing an in- going wave and the other describing an 
out-going wave. As is well-known, to obtain the retarded Green's function, we must use an 
in-going wave. We still need to choose the normalization Ua+u- As shown in refs. ||9|, [Tl[| , 
for fermions, once we choose an in-going wave solution, if we use the projectors ni^2 5 then 
when we fix the normalization of the up component ^aiu to be Uaiu, the equation of motion 
fixes the down component ^aid to have normalization i times Uaiu- The same applies to the 
up and down components of "^a2- Switching to the II-i- projectors (recall eq. ( 2.13D ), the 



statement is that once we fix the normalization of ^a-d to be ria-d, then ^a+d must have 
normalization i times Ua-d- The same statement applies to ^a-u and '^a+u- 

We thus need only fix 2A^ normalizations, for the up and down components of the ^a-- 
Let us arrange these normalizations into a row vector n 

n = {ni-u, rii-d, n2-u, n2-d, ■■■ , riN-u, riN-d) ■ (2.34) 

Following refs. [^], we use these horizon normalizations to construct a basis of solutions 
as follows. We solve the equations of motion 2A^ times, each time with a different choice 
of n. The first time we use n = (+1, +1, +1, . . . , +1, +1), the second time we use n = 
(+1, —1, +1, . . . , +1, +1), the third time we use n = (+1, +1, —1, . . . , +1, +1), and so on. We 
label these choices n^^\ with i = 1,...,2N. For each choice of normalizations, we obtain 

(i) 

solutions ^a±- have a basis of solutions, so we can write any particular solution as 

a linear combination of these. To do so, we construct matrices that we will call P^{u,uj, k) 
from the basis solutions, where each row corresponds to a field and each column corresponds 
to a choice of normalization (the i index). For example, (suppressing the ^a-'s dependence 
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on all variables) 



(2.35) 



with P^j{u,uj,k) defined similarly. The P^{u,uj,k) are 2A^ x 2N matrices. For later con- 
venience, we will factor out the leading asymptotic behavior of the solutions, defining new 
matrices P^-{u,uj,k), 

P± (n, a;. A;) ee ^.f ±- P± (^z, a;, fe). (2.36) 
We can now write any solution as a linear combination of the basis solutions: 

^a+{u,u:,k) = ui+^ P+{u,u:,k) {P+i€,u,ky')^^ Cb+iu;,k), 

^a-{u,uj,k) = J-"'p-jiu,u;,k) (p-(e,a;,/c)-i)., Cb„(u;,fe), (2.37) 

with a summation over the j index. Notice that we take the solutions ^a± to be linear in the 
"sources," Ca±- As emphasized in refs. eq. (|2.37| ) is simply saying that the sources Ca± 
will source various linear combinations of fields in the bulk, and that we can write those linear 
combinations as linear combinations of our basis solutions. Notice that when we evaluate the 
solutions at u = e, we reproduce the leading asymptotic form, ^a± ~ Ca± U2^^. 

Now we arrive at the main difference between bulk fermions and bulk bosons: Ca+ and 
Ca- are not independent. The equation of motion relates them |62, |6^. Indeed, we saw above 
that only the Ca- are sources, while the Ca+ give one-point functions (roughly speaking). To 
relate them, we follow refs. |10, We return to the equation of motion as written in eq. 



( |2.32 ). We focus only on the equation with on the left-hand-side, and simply insert 
solutions as written in eq. ( 2.37| ) (suppressing all u! and k dependence) 



V 



ab+U2+"^P+iu) 



M, 



ae+ 



U'2 



p- 



(2.38) 



where the parentheses separate u-dependent factors from ^-independent factors. We now 
observe that the matrices P^- also solve the equation of motion, by construction, since they 
are built from solutions. We thus have 



V.,+ ^^+'"P+(n) 



Mac+U^-'P-iu) 



(2.39) 



Here we have a free j index, so we actually have 2N such equations. (Recall that the index 
j labels the choice of normalization vector n.) The above equation is just the statement that 
one column of the matrices solves the equation of motion. We are free to act on the right 



with the vector (P"'"(e) ^) Crf+, so that we obtain 



V„,+ ^f+-P+(n)((P+(6)-l)^., C,+ 



Ma 



-^P-^{U) ({P^{e)-^).^ c,A . (2.40) 
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We now simply compare eqs. ( p.38|) and ( 2.40 ). The left-hand sides are identical, so we may 
equate the right-hand sides. Acting on the left with some inverse matrices, we obtain the 
desired relation between the Ca+ and Ca~, 



la] 



ijb 



Cb-. 



(2.41) 



Invoking eq. ( ^.281) , we now just need to perform two operations to extract the retarded 



P 



Ijb- 



we take e ^ and then act on the 



two-point function G^f^{uj,k) from P^{e)aj 
right with —ij^. 

The effect of taking e — )■ is easy to understand. From the definition of the P^{u) in eq. 



( |2.35| ) and the definition of the P^{u) in eq. ( |2.36| ), we can identify the e — t- limit of the 

P±(n) as 



lim 



.(1) J2) 

.(1) J2) 

-2- ^2- 



■N- 



.(2) 



"N- 



(2.42) 



and similarly for lime-s-o ^aj(e)- In short, the matrices P^-, when evaluated at the boundary. 



are simply matrices of the Ca+ and c^-. 

Notice that the P~(e)~^ matrix will introduce a factor of detP~(e) in the denominator of 
the Green's function. Generically, then, if detP~(e) has a zero, the Green's function will have 
a pole, which means a quasi-normal mode appears in the bulk spectrum, as in the bosonic 
cases of ref. ||55|| . Given the identification in eq. ( p. 42]) , then, to identify quasi-normal modes 
we need only identify the zeroes of the matrix of c_ 's. 

Understanding how 7* acts on P'^{e)aj (P~(e)~^)^.^ is a little tricky. Luckily, the way we 
have written P^{e)aj (P^(e)^^)_^.^ means that —ij^ acts trivially. To see that, notice that eq. 
( |2.41| ) is written in a two-component form: here Ca± are two component spinors. To restore 
them to four-component form, we take a direct product. 





( \ 









Ca-d 

\ y 



(2.43) 



To restore the P^^ matrices to the same four-component form, we recall recall eq. ( 2.35D , which 



shows that we should perform exactly the same direct products (suppressing the dependence 
on all variables): 

/ o\ „_ „_ / 1\ 



p: 



1 



Pal^ 



0, 



(2.44) 



which implies {P 



(P- 



/]a 



1 ) . Eq. ( ^ ) thus becomes 




P+{e)a] {P-{e)-') 



jb 





(2.45) 
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We now simply observe that, in such a representation, — ^7* = Iat^cti. In the NxN subspace 
we want, — Z7* merely acts as the identity. 

In summary, the retarded Green's function for coupled bulk fermions is 



G« (u;, k) = lim (p+(6),, P-{e)-,') 



(2.46) 



with the matrices defined in eq. ( p. 36 ) . 

Finally, as an important check, let us use our prescription to reproduce the result for free 
fermions, eq. ( p. 30 ). For illustration, we consider = 2, so we have two bulk fermions, which 
we will call and ^f,. We return to the equation of motion as written in eq. ( p. 32 ), and 
assume the equations for and ^'^ decouple, so that Vab± and Mab± become diagonal in the 
a and b indices. We can further decouple the equations of motion by using the projectors ni^2- 
Acting with these, we obtain equations similar to eq. (2.14). We thus find four decoupled 
equations, for ^ai, ^a2, ^bi and ^'62- 

We now solve the equations 2N = 4 times, each time with a different normalization 
vector n for the ^a- and ^b- fields. In the first solution, all four fields have normal- 

(+1, +1, +1, +1). In the second solution, we use 



izations n 



[n-a—ui ^a—di n-b—ui n-b—d) 



n = (+1, — 1, +1, +1). The key observation is that the field ^a-d whose normalization we 
change is ^a-d = ^aiu (recall eq. ( 2.13D ), and hence couples only to ^aid = ^a+d- The 
change in normalization thus leaves the other three fields, ^a~m ^b-m and ^b-d unchanged. 
The solutions for these fields will thus be identical to what they were using the original +1 
normalizations. The P~- matrix thus takes the form (here we must write the up and down 
components explicitly) 



f^\'^ vl;« ^(1) \ 

1— U 1—U 1—U 1—U 

$(1) ^(2) (1) (1) 

^1-d ^1-d ^1-d ^1-d 

Vl/W ,1/(3) M/W 

2— u 2—u 2—u 2—u 

\^2-d ^2-d ^2-d ^2-d, 



(2.47) 



with P^j{u,uj, k) being identical except all — subscripts become +. The main feature here 



is that all the superscripts are the same, except on the diagonal. A straightforward exercise 
(especially simple for 2x2 matrices) then shows that taking the inverse {P~)ja and then 
contracting with Pj^, and taking e — t- 0, reproduces exactly the purely diagonal Cj^u/c-u and 
c+d/c^d form of eq. ( |2.30| ). 

In summary: by combining the methods of refs. 55] and 



pq ], we have provided 

a relatively simple prescription to compute the matrix-valued retarded two-point function 
from bulk solutions for coupled fermions. We simply solve the equations of motion (typically 
numerically) 2A^ times, using a different normalization vector n each time, use those solutions 
to construct the matrices P^(e), and then take Ymii^^o P'^ {e)aj (P~(e)~"^) . 
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3. Probe Branes and Holographic p-wave Superfiuids 



In this section we review how to obtain a holographic p-wave phase transitions from simple 



string theory constructions of intersecting Dq-branes and Dp-branes 33, 35 1 (see also 
ref. |^6|). We also present some new results for the particular D3/D5 system we subsequently 
explore in later sections. 

3.1 p-waves, Probe Branes, and Vector Meson Condensation 

The minimal ingredients for a holographic p-wave phase transition are gravity in a black hole 
spacetime with holographic variable u (and some dual field theory), plus non-Abelian bulk 
gauge fields. We will consider the simple example of SU (2) gauge fields ^4^^, where a = 1, 2, 3 
labels the generators Ta of SU (2), although other non-Abelian groups besides SU{2) work just 
as well [|^, As described in the introduction, the p-wave superfluid transition appears in 
the bulk as a charged black hole growing vector hair at low temperature. 



As observed in refs. |3^, |3j, we can easily obtain holographic p-wave superfiuids 
using well-known intersections of Nc coincident Dq-branes with Nj coincident Dp-branes in 
type II string theory, which we will refer to as Dq/Dp systems. The idea is to take the 
usual decoupling limit for the Dq-branes, which in particular means Nc ^ oo, to obtain type 
II supergravity in the near-horizon geometry of the Dq-branes. Starting with non-extremal 
Dq-branes produces a black hole spacetime. 

If we keep Nj fixed as A'^c — ^ OO; so that Nj <ti Nc, then we may neglect the effect of 
the Dp-branes on the supergravity fields. The Dp-branes are then probes of the background 
geometry, and their dynamics is described by the non-Abelian Born-Infeld action (possibly 
plus Wess-Zumino terms) with gauge group U{Nj). If we introduce exactly Nf = 2 probe 
branes, then we have U{2) gauge fields on the worldvolume of the Dp-branes. The SU{2) 
subgroup gives us the SU{2) gauge fields we want. Given that the non-Abelian Born-Infeld 
action is not known to all orders in the field strength, we are typically limited to working at 
the leading non-trivial order, which is the Yang-Mills term. We have thus obtained SU{2) 
gauge fields in a black hole spacetime. 

Crucially, notice that these Dq/Dp systems give rise to probe gauge fields, rather than 
gauge fields coming from the supergravity sector. The probe limit is sufficient to study many 



properties of the p-wave phase transition pq, 32, 33 , 34, p5[, however, the probe limit is known 



to fail at low temperatures, because the solutions with nonzero Al,(u) have a field strength 
that increases as we cool the system, so that we can no longer neglect the back-reaction (or 
trust the Yang-Mills approximation to the non-Abelian Born-Infeld action) [ p2| , 34 1 . To reach 
zero temperature, we must solve the fully coupled equations of motion, which to date has 



only been done in ad hoc models [17, [l^, 23 



The benefit of the Dq/Dp construction is that we can identify the dual theory, which 
we will call the Dq/Dp theory. The Nc Dq-branes generically give rise to an SU{Nc) gauge 
theory with fields only in the adjoint representation of SU{Nc). Open strings from the Nf 
Dp-branes to the Dq-branes give rise to fields in the fundamental representation of SU{Nc), 
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i.e. flavor fields. In analogy with (supersymmetric) QCD, we will call any flavor fermions or 
scalars "quarks" or "squarks," respectively. If the Dp-branes do not overlap with all q spatial 
directions of the Dq-branes, then the flavor fields will be confined to propagate along some 
defect of nonzero codimension. 

In the field theory, the probe limit consists of neglecting quantum effects due to the flavors, 
such as the effect on the running of the coupling, because such effects are parametrically 
suppressed by Nj/Nc- In the language of perturbation theory, the probe limit consists of 
discarding all diagrams involving quark or squark loops. 

If we separate the D-branes in an overall transverse direction, we may give the Dq-Dp 
strings a finite length and hence the flavor fields a finite mass, although in this paper we 
consider only massless flavor fields (unless stated otherwise). 

The U{Nf) gauge invariance on the Dp-branes' worldvolume is dual to a U{Nf) flavor 
symmetry, analogous to the vector symmetry of QCD. The overall U{1) we identify as baryon 
(or really quark) number, and the SU{2) subgroup we identify as isospin. 

We can thus easily see what the bulk transition looks like in the field theory. We have a 
strongly-coupled, large- A^c non-Abelian gauge theory coupled to Nf = 2 species of massless 
flavor fields, which may be confined to a defect. We study thermal equilibrium states with 
temperature T, and introduce an isospin chemical potential /i for C/(l)3. For sufficiently large 
/i, the system develops a nonzero (Jf). The operator Jf is a gauge-invariant bilinear in the 
flavor flelds, valued in the adjoint of SU{Nf). For example, in two-flavor massless QCD, with 
up and down quarks u and d, Jf ~ u^^d. 

Such an operator is precisely what we would call a vector meson, and the phase transi- 
tion appears to be vector meson condensation. To be precise, the spectrum of the Dq/Dp 
theory includes gauge-invariant bound states of flavor flelds. We will refer to such bosonic or 
fermionic bound states as "mesons" or "mesinos," respectively. For massive flavor fields, these 
mesons/mesinos are typically the lightest fiavor degrees of freedom in the theory [6^, 7C]. We 
may thus imagine writing an effective theory for these degrees of freedom, analogous to the 
chiral Lagrangian of QCD. An isospin chemical potential /i will act as a negative mass-squared 
for any mesons/mesinos charged under [7(1)3. If we make fi sufficiently large, then we expect 
Bose-Einstein condensation of mesons. In QCD, we expect the lightest charged mesons, the 
pions, to condense first producing a scalar condensate (and hence an s-wave superfiuid), while 
the heavier vector mesons may condense at higher fi |7^, |7^. Which mesons condense first 
in Dq/Dp systems depends on the details of the system. In Dq/Dp holographic models of 
QCD, such as the Sakai-Sugimoto model [73|, holographic calculations suggest that indeed 
the pions condense first and the vector mesons second, as we increase ^ [74, |7^. The general 
lesson from these Dq/Dp systems is that the p-wave superfluid phase transition appears to 
be vector meson condensation, which is in line with our weak-couling intuition. 

Moreover, thinking of the p-wave states as a Bose-Einstein condensate makes many poten- 
tially confusing features of the p-wave state transparent. For example, the p-wave transition 
appears to involve the spontaneous generation of a persistent current (Jf ), that is, at high 



- 20 - 



density charges begin moving without experiencing dissipation.^^ While not impossible, such 
a scenario naturally raises some questions. Why do charges start moving? How does that 
lower the free energy? Vector meson condensation neatly accounts for all of the physics: we 
merely see Bose-Einstein condensation, i.e. bosons populating a zero-momentum state, the 
main novelty being that the bosons are vectors, not scalars. 

We are interested in condensed matter applications, and in particular quantum critical 
theories, which are scale-invariant, hence our Dq-branes will be D3-branes, whose near- horizon 
geometry is ^^5*5 x S^. The dual theory is then a CFT, namely (3+l)-dimensional A/" = 4 
SYM with large Nc and and large 't Hooft coupling. Our Dp-branes will preserve half the 
supersymmetry of the background, which means Dp-branes extended along AdSp x with 



P + Q = p + 1, where supersymmetry requires |P — Ql = 2 |36]. Well-known examples include 



D7-branes extended along AdS^ x 5^ ||7|, ^ or AdS^ x S'^ @ |7|, D5-branes extended 
along AdS4^ x 5^ |2|, or D3-branes along AdSs x fH-^^ For two coincident D7-branes 
with P = 5 and Q = 3, holographic calculations have shown that a p-wave transition occurs 
precisely when peaks in the Green's function, namely those corresponding to vector mesons 
charged under U{1)3 in the vacuum state, cross into the upper-half of the complex frequency 
plane, indicating an instability toward Bose-Einstein condensation. 

With an eye toward condensed matter applications, and for technical reasons we explain 
in section ^, we will work with (two coincident) D5-branes with P = 4 and Q = 2. The 
dual theory is thus = 4 SYM coupled to Nj = 2 flavor fields that propagate only in 2+1 
dimensions (a codimension one defect) and preserve (2+l)-dimensional = 4 supersymme- 
try (eight real supercharges). The field content and Lagrangian of the D3/D5 theory were 



determined in refs. [|4l| , |42|| , which we review in section The D3/D5 system exhibits rich 
thermodynamics, studied in detail (holographically) in |43, |4^, |4^, 47, |5^, 51]. Here 
we initiate the study of the D3/D5 theory with a finite isospin chemical potential and finite 
temperature. As expected, we will find a p-wave phase transition. 

Our ultimate goal is to compute, holographically, fermionic retarded Green's functions 
in p-wave superfluid states. We discuss fermionic excitations on the worldvolume of the 
D5-branes, and the dual mesinos, in section ^ We will be able to compare our numerical 
results with previous studies, however, we cannot exploit the analytic results for the form of 
fermionic Green's functions derived in ref. The analysis of ref. involved fermions 

in an extremal Reissner-Nordstrom background, and in particular made great use of the 
emergent near-horizon AdS2 factor. Without having access to the T = finite-density state, 



Crucially, however, no net momentum is flowing. In holographic calculations, in both the probe and fully 
back-reacted cases 0, ^ H, the Yang- Mills stress -energy tensor and the metric are diagonal, and indeed 
the bulk spacetime is static, which indicates that the expectation value of the field theory stress-energy tensor 
is strictly diagonal. The system thus has zero net momentum. If charges are moving, they must be doing so 
in pairs that move in opposite directions. A static bulk spacetime also indicates that the energy density of the 
field theory is not changing in time: the system is not heating up, consistent with the fact that the moving 
charges experience no dissipation (no frictional forces). 

^^These are the cases in which the flavor fields propagate in at least one spatial dimension, and hence a 
nonzero (Jf ) is possible. 
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we do not know for sure whether the geometry exhibits an emergent AdS2 or something 
similar. Nevertheless, the equations of motion for our fermions are formally similar to those 
of refs. |ll|, and hence we can recover similar finite-temperature results. 

Despite the bad news that we cannot reach T = in the probe limit, we do have good 
news: we can study fermionic response near the p-wave transition. The high-temperature 
normal phase is rotationally symmetric, but the p-wave phase is of course not. We will see 
the breaking occur explicitly in fermionic spectral functions in section ^ 

3.2 Probe Dp-branes in AdS^ x 

We now want to study a p-wave superfluid transition for (2-|-l)-dimensional flavor fields 
described holographically by two coincident D5-branes with P = 4 and Q = 2. Although our 
(numerical) analysis will be for the D5-brane, in the interest of generality, and to connect to 
our discussion in section |^, we will write formulas for an arbitrary Dp-brane with P > 3. 

The background supergravity solution includes a metric and Ramond-Ramond (RR) five- 
form. The five- form will be important in section ^ The spacetime is (4-1-1 )-dimensional 
AdS-Schwarzschild times S^, with metric 

with 



f{u)dt'^ + dx^^ +dsl5, (3.1) 



f{u) = 1 - 74 , nh = —. (3.2) 



^ -1. 
< 

In our units, where the AdS radius is one, we can convert between string theory and field 
theory quantities using q'~^ = AttQsNc = 2gY]^jNc = 2A. 

As we will be studying fermions in section^ we will need the vielbeins and spin connection 
associated with the metric above. We record these here for later use. The nonzero vielbeins 
of (4-|-l)-dimensional AdS-Schwarzschild are (recall that upper index is local Lorentz, and 
the lower index is general coordinate), 

e\ = -^, e\ = ^, e'j = -S'^. (3.3) 

The spin connection uj of (4-|-l)-dimensional AdS-Schwarzschild then has the nonzero com- 
ponents 

uJtu = — dt, = ax, (3.4) 

\ u 2 J u 

where uj^u indicates the three components oJxu-, ^yu, and ujzu- 

Next we introduce two coincident probe Dp-branes extended along AdSp x S^. We will 
only consider the trivial embedding of the Dp-branes, that is, we consider solutions in which 
all the Dp-branes' worldvolume scalars (including scalars in AdS^ directions) are zero. The 
dual flavor fields are then massless. The induced metric on the Dp-branes 
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where now dx'^ represents the appropriate-dimensional EucHdean metric. 

We want non-trivial worldvolume SU{2) gauge fields. The action for the gauge fields, to 
leading non-trivial order, is 



Sdp = -TopNf J cp+^i 



l + (2W)2iTr(F^,F^-) 



where T^p = {2tt)~^ (a') ^ is the tension of the Dp-brane, Nj = 2, the integral is over 
the worldvolume coordinates gop is the determinant of the induced metric, and the trace 
is over gauge indices. We use SU{2) generators Ta = such that, with €123 = +1 

[Ta,n] = ieabcTc- (3.6) 

The field strength Ff^^, = F^^ Ta- 

The equation of motion for the gauge field is simply the Yang-Mills equation, 

v^Fr+/afec(^^)6i^r = o. (3.7) 

For probe Dp-branes wrapping AdSp C AdS^ with P > 3, we will consider solutions of eq. 
( |3.7D of the form 

A = Al{u)Tidx + Al{u)T2.dt, (3.8) 
in which case the Yang-Mills equation becomes 

{4)" + (A^)' - -j^-A^ {Alf = 0, {3^9a) 

(-^i)" + (^ + 7f^) (-^9' + IW- = 

where primes denote du- The equations of motion determine the asymptotic forms of the 
solutions, 

Al{u) = fi-d^ u^-^ + . . . , Al{u) = dl u^-^ + . . . , (3.10) 

where . . . represent terms that decay faster than u^~^ as ti — t- 0. Here the constant d^ is 
related to (J3) as 

(J*) = NfTop {2^a'f df = {2tt)'P+^ 2^ NfN^X^ d^, (3.11) 

where in the second equality we converted to field theory quantities. Similarly, (Jf ) = 
NjTdp {2TTa')'^ d^.. Notice in particular that both (J|) and (Jf) are proportional to NjNc- 
One solution of eqs. (|3.9D has Al.{u) = and 

AUu)=fi\l--p^j . (3.12) 



^^Notice that Al.{u) has no leading constant the way j4|(m) does, so no source for Jf is present in the field 
theory: (7(1)3 will be broken spontaneously. 
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Such a solution corresponds in the field theory to the normal phase, in which the chemical po- 
tential n explicitly breaks the SU{2) isospin symmetry down to the U (1)3, but no spontaneous 
symmetry breaking occurs. These solutions exist for all values of ^. 

Notice that these D3/Dp theories in the probe limit are scale- invariant^^, so the only 
meaningful physical quantity is n/T, so fixing T and increasing fj, is equivalent to fixing 
/i and reducing T. We will think in terms of the latter. Any transition must occur at a 
temperature Tc set by the chemical potential, Tc oc /i. 

For sufficiently low T (or large fi), other solutions of eq. ( p.7[ ) exists in which Al.{u) is 
nonzero. For the D7-brane with P = 5 and Q = 3, such solutions were found numerically 
in refs. These solutions correspond in the field theory to superfluid states, 

with nonzero (J^), so ^7(1)3 is spontaneously broken. For P > 4, the field theory's spatial 
rotational symmetry^^ is also broken from SO{P — 2) down to SO{P — 3). 

For sufficiently low T, we have two solutions, so we need to determine which is thermo- 



dynamically preferred. As shown in refs. p2| , 33, 34| for the D7-brane with P = 5 and Q = 3, 



the superfluid phase is thermodynamically preferred relative to the normal phase for all T/ii 
where the solutions with nonzero Af{u) exist. The transition between the phases is second 
order, with mean-field exponents. In particular, near the transition, the condensate has mean 
field exponent 1/2: (Jf ) oc (1 - T/T^f^'^. 

For two coincident D5 branes with P = 4 and Q = 2, the story is qualitatively the same. 
For fj. > 3.81 X (vrT) the state with nonzero (J^) has lower free energy. In other words, 
Tc = 3 g^xTT - figure we plot the constant dl., which is proportional to (Jf), versus the 
rescaled temperature T/Tc- Near the transition, (Jf ) appears to have a mean-field exponent 



of 1/2 as in the D7-brane case [l32| , p3| , 34 1. 



Notice that as T decreases, the condensate grows, which, as explained in refs. |14, 15], 
suggests that we are leaving the probe limit. We will only present results in the p-wave 
phase for T/Tc > 0.4 (as plotted in figure ([l|)), where we have some hope that the probe 
approximation captures the essential physics faithfully. 

4. The World volume Fermions 

We want to study fluctuations of fermionic operators of the D3/D5 theory with finite tem- 
perature and isospin chemical potential, in the two phases described above, the normal (non- 
superfluid) phase and the superfluid phase. On the field theory side, we will study mesino 
operators valued in the adjoint representation of the SU{2) isospin symmetry. As we will see 
below, we will thus have three mesinos, two with equal and opposite charges under the U{l)s, 



and one that is neutral. We discuss the form of the mesino operators in section 4.2 



^^In the probe limit we neglect the quantum effects that would cause the N' — 4 SYM coupling to run, and 
we are working with massless flavor fields. The theory is thus in a limit where no intrinsic scale appears. 

^^When P = 3 the flavor fields are conflned to a (H-l)-dimensional defect and hence have no spatial 
rotational symmetry. Notice that in the P = 3, 4 cases the large Nc limit is what permits spontaneous 
symmetry breaking to occur, by suppressing the fluctuations that would destroy long-range order. 
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Figure 1: The condensate (J^), times jf-j^^j^, versus the rescaled temperature T/T^- 



To be specific, we will compute holographically the retarded two-point function of mesinos 
as a function of frequency and momentum. On the gravity side, that means studying fermionic 
fluctuations of the Dp-branes, and in particular solving their linearized equation of motion 
in the background we found in the last section, where the geometry is (4+l)-dimensional 
AdS-Schwarzschild and the D5-branes have non-trivial worldvolume gauge fields. As for 
all super symmetric Dp-branes, the worldvolume fermions are in a supermultiplet with the 
worldvolume gauge field and scalars, and hence are in the adjoint of the worldvolume SU{2) 
gauge group, which is dual to the statement that the mesinos are in the adjoint of the isospin 
symmetry. 

As we will see, in the normal phase the three bulk fermions decouple, which is dual to the 
statement that the retarded Green's function is diagonal. In the superfluid phase, however, 
where Al.(u) is nonzero, the three fermions couple, indicating that the dual operators mix 
under renormalization group flow in the field theory. The Green's function is then a 6 x 6 
matrix, where the 6 is the number of fermions times the two components of the fermions 
(two, using the ni^2 projectors). We thus have a perfect testing ground for the method we 
developed in section |2| for computing Green's functions from coupled bulk fermions. 

The fermionic part of general Dp-brane actions, to quadratic order in the fermionic fields 



and in backgrounds with non-trivial RR forms, was determined in refs. 38, The 

general couplings were derived by starting with the action for a supermembrane in M-theory, 
written in a superspace formalism, expanding the action to second order in the Grassmann 
variables, reducing to type IIA supergravity, and then performing a T-duality to type IIB. 



^^Obviously, the mesinos carry no baryon number, which is dual to the statement that the worldvolume 
fermions, like all of the worldvolume fields, do not couple via a gauge-covariant derivative to the diagonal (7(1) 
part of the worldvolume gauge field. 



The form of the quadratic fermionic action on the Dp-brane worldvolume is thus determined 
completely by supersymmetry and T-duality. 

Using the worldvolume fermion actions of refs. |3^, the spectra of mesinos in the 

D3/D7 theory (for the D7-brane with P = b and (5 = 3) and in the Sakai-Sugimoto model 
were determined in refs. |4^, We will very closely follow the D7-brane analysis of ref. 
[^ ] , which in turn was the fermionic generalization of the analysis of ref. |6^] for mesons. For 
a Dp-brane extended along AdSp x , we consider a worldvolume spinor that is a spherical 
harmonic on the S"^. We reduce the worldvolume spinor on the 5*^, obtaining an effective 
Dirac action in AdSp. This procedure fixes the masses of the bulk fermionic excitations, 
which allows us to identify the dimensions of the dual mesinos, and more generally to map 
bulk fluctuations to mesino operators. As emphasized in ref. |4^, the coupling to the RR 
five-form is crucial to obtain the correct bulk masses. 

One of our main points is: because we work with a particular string theory system, we do 
not have the freedom to change the mass or the charge of our bulk fermion, in stark contrast 



to the models of refs. P, IC, O]. Both the masses and charges are, ultimately, fixed by 



supersymmetry and T-duality, as explained above. 
4.1 Equation of Motion I: Reduction to AdS 

We will now repeat the analysis of ref. |40|, in which the fermionic action of a D7-brane 
extended along ^^5*5 x was reduced to an effective Dirac action in AdS^,, but now for more 
general Dp-branes extended along AdSp x 5"^, with emphasis on D5-branes with P = 4 and 
Q = 2.21 

The quadratic action for fermionic fluctuations of the Dp-branes is (refs. ||3^, ^) 

Here ^' is a ten-dimensional positive-chirality Major ana- Weyl spinor of type IIB supergravity, 
the are the pullback of the ten-dimensional T-matrices to the Dp-brane worldvolume, 
^ A ~ ^ M^A^^ ^ trivial embedding, so the pullback is trivial), P_ is a K-symmetry 

projector that ensures /t-symmetry invariance of the action. Da is a (gauge and curved- 
space) covariant derivative, and F^p^j^^ is the five-form of the background. Notice that 
here Fjqpqfis is not the pullback of the five- form to the Dp-brane worldvolume, rather, it 
is the five-form evaluated on the submanifold spanned by the Dp-brane. Indeed, no part of 
the expression FfjpQps^^^^^^ involves a pullback. Here A,B,... denote all worldvolume 
indices, while below A^B,. . . denote AdS-Schwarzschild coordinates which are wrapped by 



^"The authors of ref. appear to omit the couphng to the five-form when they study fermionic fluctuations 
of the probe D3-brane extended along AdS^ x . 

^^In section ^ we were interested in p-wave states and hence required P > 3. In this subsection we relax 
that constraint. Our results will thus also apply for D5-branes extended along AdS2 x 5"*, which were used in 
ref. 184] to construct a holographic model of fermions at lattice sites (that can pair to form dimers). 
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the probe brane. Moreover, the indices of the coordinates on the sphere S'^ are labelled by 
a,b, . . . . Notice that the fermion ^ is in the adjoint representation of SU{2). 

The equation of motion for the fermion is (for now we suppress gauge indices) 



^ 8 2*5! NPQRS y 



^ = (4.1) 



We will reduce the equation of motion for the fermion to a Dirac equation in AdSp, following 
ref. [ pO| ] very closely. First we decompose every ten-dimensional spinor and F-matrix into 
parts associated with AdS^ and S^. In a local Lorentz frame, the F-matrices decompose as 

F^ = ^ I4 ^ 7*^, F™ = ai0 7™®l4, (4.2) 

where I4 is the 4x4 identity matrix, the index M runs over AdS^ directions (which we will 
generically call 01234), and the index m runs over directions (which we will call 56789). 
The 7-matrices are five-dimensional, obeying the usual relations 

{/^,7^} = 2r/^^^, {7^,7"} = 25"^". (4.3) 

Given the above decompositions, we then have 

F01234 ^ ^14^ I4, r^'^^S^ = CTi ® I4 ® I4 (4.4) 

pll^p0123456789 ^^^g^l^^-^^ (4.5) 

The ten-dimensional spinor ^ has positive chirality, F^^^ = ^, and decomposes as 

^'=t«>X«>^, (4.6) 

where t= ^ q ^ ' ^^'^ ^ ^^"^ ^ four-component spinors of SO{5) and S'0(4, 1), which act 

on the tangent spaces of and AdS^, respectively. The spinor x further decomposes as 
X = X\\ ® X±) where x\\ is a spinor associated with the S*^ that the Dp-brane wraps and x± 
is associated with the directions transverse to the S^. 

We parameterize the five- form in terms of the volume forms of AdS^ and S^, which we 
denote as i^AdS^ and 



FnPQRS — 4 i^AdSs) NPQRS ' -^npgrs — 4 (ri^s. 



npqrs 



Using the decomposition of the type IIB spinor in eq. (4.6) and of the 7-matrices in eq. ( [4. 21 ) 
we obtain 

pa (p01234 ^ p56789^ F„ (f ®X ® ^) = 2 i ®X ® ^ , 
pA ^p01234 ^ p56789^ F^ (t ®X ® ^) = "2 i ®X ® ^ , (4-7) 



-27- 



In eq. (4.7) we have not summed over a or A. Using this result we can simphfy the couphng 
of the spinor to the five-form, 

\i^^^Fnpqrs^^^^^''^a^ = ((^1 + ^^2) ®U® I4) r^^- = -\{P-Q) (i ®x ® ^) 

(4.8) 

where here we do sum over A. We can also extract the S'^ and the AdSp part of the derivative 
terms as 

= {{<T2®U® i'^Da) + (0-1 ® -i'^Da ® I4)) (t ®X ® ^) 

= {i (I2 ®U® i'^Da) + (I2 ® 7" Da (g) I4)) (i (E)X ^) 

^ {iPAdSp + PsQ)ii^X^'^), (4.9) 

where PAdSp ^'^d PgQ are the Dirac operators of AdSp and S'^ , respectively. The Dirac 
operator on a sphere has spinor spherical harmonics that obey 

+ _ ^ ^ ( n , Q\ + 



^sQxt = ^^\^^+^jxt, (4.10) 

where ^ > and Rq is the radius of the . In our units, Rq = 1. For Q = 3, relevant for 
the D7-brane along AdS^ x S^, the spinors xl are in the (^^, |) representation of the 50(4) 
that acts on , while the spinors x7 ™ the (|, ^^) representation. For (5 = 2, relevant 
for the D5-brane along Ad5'4 x 5^, the spinors are in the + ^) of the SU{2) ~ 50(3) 
that acts on 5^. 

Inserting everything into eq. (|4.1| ) , we find 



The fermions^^ thus have masses (in our units, where the radius of AdS is one) 

m+ = l + ^, mj = -{^ + Q-]^P^. (4.12) 

We collect the values of for our Dp-branes of interest in the table below. 

Notice that since P and Q are integers, the rri^ will always be integer or half-integer. As 
we review in the appendix, a bulk fermion with an integer or half-integer mass m is dual to a 
fermionic operator of dimension A = + |m|. We include the values of = + \ni^\ 
in the table. 

To get a rough idea of which operators correspond to which bulk fermion, we can do 
some dimension counting. Let us denote a generic quark as V') a generic squark as a generic 



^■^The ± superscript here refers to the sign of the eigenvalue of the Dirac operator in eq. (4.1C), not to the 
projectors n± defined in section ^. 
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Dp 


P 


Q 


m 




mj =i + P/2 


A/ 


1^71 =£ + Q-P/2 




D7 


5 


3 


3/2 


1 


£ + 5/2 


£ + 9/2 


£+1/2 


£ + 5/2 


D7 


3 


5 


1/2 




£ + 3/2 


£ + 5/2 


£ + 7/2 


£ + 9/2 


D5 


4 


2 


1 


1/2 


£ + 2 


£ + 7/2 


£ 


£ + 3/2 


D5 


2 


4 







£+1 


£ + 3/2 


£ + 3 


£ + 7/2 


D3 


3 


1 


1/2 





£ + 3/2 


£ + 5/2 


£-1/2 


£+1/2 



Table 1: Masses of fermionic excitations on the worldvolume of a Dp-brane extended along AdSp x S'^ 
inside AdS^ x . We list Dp-branes that are known to preserve eight real supercharges (at zero 
temperature and density), in which case |P — Q| = 2. Here if^ denotes a generic quark field and q 
denotes a generic squark field. denotes the dimension of the operator dual to the bulk fermion 
with mass m^, with = -^^^ + \mf\. For the D3-brane, the values of mj and A^ shown are for 
£ > 1 only, whereas for £ = 0, Img | = 1/2 and Ag = 3/2. 



adjoint Majorana fermion as A, and a generic adjoint real scalar as X. The dimensions of the 



fields are [■0] = [q] = [A] = 3/2, [X] = 1. 



For the D7-brane extended along AdS^ x S^, the D5-brane along AdS^ x 5*^, and the 
D3-brane along AdS^ x S^, all of which have P — Q = 2, the dual flavor fields comprise 
a supermultiplet with both quarks Tp and squarks q. In these cases, we can build a gauge- 



invariant mesino in two ways [40|. One way is to construct an operator of the form t/^A^' + 
q^XXq, with dimension A = P — 1/2. We can additionally include some number £ of adjoint 
scalars^^ as tpXX^il) + q^XXX^q, so that the dimension is A = £ + P — 1/2. Inspecting the 
table, these are precisely the A^, so apparently these kinds of mesinos are dual to the bulk 
fermions with masses m^. The other way to build a mesino is to construct tpX^q (plus the 
Hermitian conjugate), with dimension A = £ + P — 5/2. For Dp-branes with P — Q = 2, these 
dimensions are precisely the Aj , so apparently mesinos of this type are dual to the fermions 
with masses mj . 

For the D7-brane extended along AdS^ x and D5-brane along AdS2 x S^, which have 



P — Q = —2, the dual flavor fields are quarks alone, with no squarks |7g, [79|, qO|, q5|. The 
mesinos with dimensions are of the same form, tpXX^il^, but the mesinos with dimensions 
Aj must obviously have a different form. We leave a detailed study of these mesinos for the 
future. 

Looking at the table, we immediately notice that the D5-branes are special: for these, 
the masses of the worldvolume fermions are integers. The reason is that the D5-branes wrap 
even-dimensional spheres, so the eigenvalue in eq. ( P^ is ±i times an integer. 

In the next subsection we focus on the D5-brane extended along AdS^ x S^, explaining 
in more detail the symmetries of the theory and the form of the mesinos. In the subsequent 
sections, we focus on the single worldvolume fermion with mass itiq = 0, for a number of 



^''Notice these are not necessarily all the same scalar, i.e. could represent £ distinct scalars. At the 
moment we are just counting dimensions, ignoring this subtlety. 
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reasons. First, of all the worldvolume fermions, these have the smallest mass, hence the 
dual operator will have the lowest dimension, Aq = 3/2, and hence be the most relevant 
mesino. Second, a numerical analysis is simpler when the fermion's mass is zero. Third, with 



a massless bulk fermion we can directly compare to the results of refs. |^, 11|, where most 
of the analysis focused on massless bulk fermions. Fourth, as we show in the appendix, a 
massless fermion requires no counterterms. 

4.2 The Dual Operators 

In this section we focus on the D5-brane along AdS4, x S'^ and study in detail the operators 
in the D3/D5 theory dual to the fermionic fluctuations considered above. For the D7-brane 
along AdS^ x S^, a similar analysis appears in ref. We begin with Nf = 1 and generalize 
to Nf > 1 at the end. 

The dual field theory is (3+l)-dimensional A/" = 4 SYM coupled to defect flavor fields 
preserving (2+l)-dimensional A/" = 4 supersymmetry (eight real supercharges). The couplings 
of the theory were determined in refs. pl| , Coupling the defect fields to the ambient fields 
requires decomposing the (3+l)-dimensional A/" = 4 multiplet into two (2+l)-dimensional 
J\f = A multiplets, a vector multiplet and a hypermultiplet. The bosonic content of the (3+1)- 
dimensional = 4 multiplet is the vector A^^ and six scalars^^ X^, , . . . , X'^ . The bosonic 
content of the (2+l)-dimensional vector multiplet is the (2+l)-dimensional vector field Ak 
and the three scalars Xy = {X'^ , X^ , X'^) . The bosonic content of the (2+l)-dimensional 
hypermultiplet is the scalar A^ and the three scalars Xh = [X/^^X^^X^). The fiavor fields 
form a (2+l)-dimensional hypermultiplet with two fermions (quarks) ifj and two complex 
scalars (squarks) q. 

The classical Lagrangian preserves (2+l)-dimensional 5*0(3,2) conformal symmetry but 
breaks the S0{&) R-symmetry down to a subgroup SU{2)h x SU (2)y , under which the scalars 
in Xh transform in the (1,0) representation and the scalars in Xy transform in the (0,1). 
We use an upper index to denote these representations: Xy and Xj^. The adjoint fermions 
A*™ transform in the (1/2, 1/2). Here i is the SU{2)v index and m is the SU{2)h index. The 
quarks transform in the (1/2, 0) and the squarks transform in the (0, 1/2). In table ^ 
(borrowed from ref. |4^), we summarize the field content and quantum numbers, including 
the conformal dimensions of the fields. 

Let us now match fiuctuations of the D5-brane probe to dual field theory operators, build- 
ing on the matching of bosonic fields and operators in refs. [^, ^]. These fiuctuations 



correspond to mesonic operators in the dual theory ||6^, 76 1, which can be arranged into a 
(2+l)-dimensional massive AA = 4 supersymmetric multiplet. 

First we consider the bosonic fiuctuations of the D5-brane, as studied in refs. |41, 7C 



The bosonic fluctuations consists of three real scalars, which in the notation of ref. [41| are 



^"^In the initial type IIB D3/D5 intersection, the D3-branes are extended along 0123, and these scalars 
represent fluctuations of the D3-branes in the 456789 directions, hence our notation. The D5-branes are 
extended along 012456, so they break the SO{6) rotational symmetry in 456789 down to SO{3) x SO{3) ~ 
SU{2) X SU{2), one rotating 456 and one rotating 789. 
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Mode 


Spin 


SU{2)h 


SU{2)v 


SU{N,) 


A 




1 








adj 


1 


V 








1 


adj 


1 


^3 











adj 


1 







1 





adj 


1 


yira 


1 

2 


1 

2 


1 

2 


adj 


3 
2 







1 

2 





N 


1 

2 


1p' 


1 

2 





1 

2 


N 


1 



Table 2: The field content of tlie D3/D5 tlieoiy. (Adapted from ref. @.) Here Ak, X^, A3, Xj^ 
and A™ are the adjoint fields of (3+l)-dimensional Af = 4 SYM decomposed into (2+l)-dimensional 
J\f = A multiplets. Ak and Xy are the bosons in a (2+l)-dimensional vector multiplet while A3 and 
Xjj are the bosons in a (2+l)-dimensional hypermultiplet. and ^/j* are the (2+l)-dimensional flavor 
fields, which are in an A/" = 4 hypermultiplet. 



Mode 


A 


SU{2)h 


SU{2)v 


Operator 


Operator in lowest multiplet 


bf 




1 + 2 


I, l>0 





Ji 


^^^^^^^^^7^^ 


(t>i 




1 + 2 


I, l>0 


1 


Si 


i^iaf-^Pj + 2q"'X^<'T<'q'^ 


{h + z)% 


-) 
-1 


l + l 


l + l, / > 





Ci 


^ I n 


(b + zt 


-) 
-1 


1 + 3 


l-l, / > 1 










Table 3: The bosonic fluctuations and their dual field theory operators for the D3/D5 system. 
(Adapted from ref. [0.) Here a are Pauli matrices, are the generators of SU{2)v, and are the 
(2+l)-dimensional T-matrices. 

4)1, (6 + z)\ ^ and (6 + z)^i^\ as well as a vector b\. Here cf^i corresponds to fluctuations of 
the embedding in 5^ directions (transverse to the S"^), {h + z)\ ^ and (6 + -z)!"*"^ are linear 
combinations of the fluctuations of the 5^ components of the worldvolume gauge field with the 
fluctuation of the embedding in AdS^ transverse to AdS^, and {k = 0,1,2) corresponds 
to fluctuations of the worldvolume gauge field in the AdS^ directions. We summarize the 
quantum numbers of these fluctuations in table ^ (borrowed from ref. ||4l|). Notice that our 
definition of I differs from that in ref. In our notation, fluctuations with the same I have 
the same mass. Later we will show that all operators with the same quantum number I fit 
into a super multiplet. Note that (6 + z)i~^^ is not present in the lowest multiplet with / = 0. 

We studied the fermionic fluctuations^^ "^f^ of the D5-brane in the last section. The 
fluctuations with I > correspond to operators with dimensions A^" = 1 + 3/2 that are 

^^In the last subsection we used a subscript £, while here we use a subscript I. For 'ifj the two are identical: 
£ = l. For *+ we take 1 = 1-1. 
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D5-brane Mode 


A 


SU{2)h 


SU{2)v 


Operator 




Z + 3/2 


1 + 1/2, l>0 


1/2 






Z + 5/2 


1-1/2, 1>1 


1/2 


Qi 



Table 4: Matching between fermionic fluctuations of the D5-brane and field theory operators. The 
explicit form of the fermionic operators J-i and Qi appear in the text below. 



in the I + 1/2 representation of SU{2)ij. The fluctuations with / > 1 correspond to 

operators with dimensions A^^-^ = / + 5/2 in the / — 1/2 representation of SU{2)h. Since 
both fluctuations are fermionic they transform in the 1/2 representation 1/2 of SU{2,)v- We 
summarize the fermionic modes in table 

We first review the lowest multiplet, i.e. I = 0, which appears already, including the 
fermionic operators, in ref. |41]. According to tables H and ^ the D5-brane fluctuation 



corresponding to the lowest-dimension operator is (6 + -z)^^^] with / = 0. Only one operator 
exists on the field theory side with the same quantum numbers as ih + z)^ \ so we can 
match (6 + z)[^ ^ with the operator Cq = q^'^aj^^q^, where are the Pauli matrices of 
SU{2)h. Co transforms in the (1,0) representation of SU(2)h x SU{2)v. Moreover Co is the 
lowest chiral primary in the multiplet since all other operators dual to D5-brane fluctuations 
have larger conformal dimensions. We can thus construct all operators in the same multplet 
as Co by applying supersymmetry generators to Cq. The supersymmetry generators form a 
2x2 matrix of Majorana spinors r/*™, which transforms like A*"*, i.e. in the representation 
(1/2, 1/2) of SU{2)h X SU{2)v. Applying the supersymmetry generators to Co we obtain the 
fermionic operator = ^/>*g™ + q^"^'tp^ with conformal dimension A = 3/2 and SU{2)h x 
SU {2)y quantum numbers (1/2, 1/2). J^q™ is dual to the fermionic D5-brane fluctuation ^^'Lq. 
Applying another supersymmetry generator to J-'q™ we obtain either j7o or £q, the forms of 
which appear in table |^. Both and £q have conformal dimension A = 2 and are singlets 
under SU{2)h but can be distinguished by their SU{2)y quantum number: Jq is a singlet 
whereas £q is a triplet under SU{2)y. 

Let us now discuss the general multiplet dual to the higher-/ fluctuations of the D5-brane. 
As in the / = case, we construct the multiplet by applying supersymmetry generators to the 
lowest chiral primary in the multiplet, C;, which is dual to (6 + z)\ \ According to ref. |41], 



the lowest chiral primary is cjo^^—^i = (jj^^o (X|^)^^"'^'^ , where (^j^) stands for the traceless 
symmetric product of / copies of the fleld Xj^. Ci has conformal dimension A = / + 1 and is 
in the (/ + 1,0) representation of SU {2)h x SU (2)y . Applying a supersymmetry generator to 
Ci we flnd the fermionic operator J^i with conformal dimension A = / + 3/2, which is dual to 
the D5-brane fluctuation . Fi is in the (/ + 1/2, 1/2) representation of SU{2)h x SU{2)y. 
Explicitly, Ti is of the form 

j,h...I,rm ^ (^X^^y^ '-^' qm + ^t™ (^X^^y'"'^' . (4.13) 

Applying another supersymmetry generator to J-i we obatin J'l or £i, which have the same 
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conformal dimension A = 1 + 2, but differ in the SU{2)h x SU{2)v representation. J'l 
transforms in the (/,0) representation whereas £i has quantum numbers 1). To obtain the 
precise form of J^i or £i we insert the operator X\j into the operator J'q or Sq, respectively. 

In contrast to the / = multiplet, other operators also appear in the multiplet for 
Z > 1, which we construct by applying three or four supersymmetry generators to C/: a 
fermionic operator Qi and a bosonic operator P/. Qi has conformal dimension A = 1 + 5/2 
and SU{2)h x SU{2)v quantum numbers (/ — 1/2,1/2). These are precisely the quantum 
numbers of the fermionic D5-brane fluctuation Explicitly, Qi has the form 

Gl =^ [^H ) X 1pj + q^ i^Xjj j A XH,l(Tnp(l ■ (4-14) 

Finally T>i has conformal dimension A = / + 3 and SU{2)h x SU{2)v quantum numbers 
(/ — 1, 0) and therefore can be identified with the D5-brane fluctuation (6 + 

We have constructed the supermultiplet for the cases / = and / > 1. For / = they 
multiplet consists of the bosonic operators Co, J^o and Sq and of the fermionic operator Tq. The 
multiplet includes are eight bosonic and eight fermionic degrees of freedom. The multiplet 
containing C;, I > 1 has 16/ + 1 fermionic and bosonic degrees of freedom: 

• One real scalar Ci in the (/ + 1, 0) representation with A = / + 1 

• One spinor J^i in the (/ + 1/2, 1/2) representation with A = / + 3/2, 

• One massive vector J'l in the (/,0) representation with A = / + 2, 

• One real scalar £i in the {1, 1) representation with A = / + 2, 

• One spinor Qi in the (/ — 1/2, 1/2) representation with A = / + 5/2, 

• One real scalar Di in the (/ — 1,0) representation with A = / + 3. 

Moreover we mapped the operators in the supermultiplet to the fluctuations of the probe 
brane summarized in tables ^ and ^ 

Finally, we consider Nf > 1 coincident probe D5-branes. The dual field theory then has 
A'^^ massless flavors, with a global U{Nf) flavor symmetry. The overall U{1) we identify as 
baryon (more accurately quark) number, while the SU (Nf) subgroup we identify as isospin. 
The mesino operators J-/ and Gi of course have zero baryon number charge and are valued in 
the adjoint of SU{Nf). For example, in our case with Nf = 2 the mesinos acquire an SU{2) 
isospin index, 77 and Q^. 

As explained at the end of the last subsection, for our numerical analysis we use the 
D5-brane fermions with zero mass, . The dual fermionic operator is J-q ^ 'ipq + q^ip. 

4.3 Equation of Motion II: Gauge Couplings 

In this section we return to the equation of motion for the worldvolume fermions, eq. ( [4.11 ), 
and specialize to our case of interest, namely two coincident Dp-branes in (4+l)-dimensional 
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AdS-Schwarzschild with trivial worldvolume scalars but non-trivial worldvolume gauge fields 
A^{u) and More specifically, we will explicitly unpack the gauge- and curved-space co- 

variant Dirac operator P^dSp the AdSp submanifold of (4-1-1 )-dimensional AdS-Schwarzschld 
and see how, when Al,{u) is nonzero, the three worldvolume fermions couple to one another. 
In this subsection we assume P > 3. 

The linearized equation of motion for the worldvolume fermions in eq. (4.11) is^^ 



AdSp 



m 



(4.15) 



where the masses appear in eq. ( 4.12 ). To simplify the notation we write m instead of 



as well as ^ instead of . PacLSp — ^ a'-^m is the gauge and curved-space covariant Dirac 
operator. The index A runs over the worldvolume directions inside AdS^. Notice that here 

are the inverse vielbeins. The 7^ are the F-matrices of (4-1-1 )-dimensional Minkowski 
space, which obey {7'^,7'^} = 2ri^^ . 

For the Dirac operator, we have (here a is a gauge index) 



AdSp 



nV77"9„ + -^7*9t + ny9^ + 
+ e^ii7^ [AM,^l-m^a, 



(4.16) 



where /' = duf- When T = and hence f{u) = 1, the operator in parentheses on the 
right-hand side in the first line is the Dirac operator of AdSp. At finite temperature, where 
f{u) = 1 — /u\, the operator is that of an AdSp submanifold of (4-|-l)-dimensional AdS- 
Schwarzschild.^^ The coupling to the gauge field in the second line is of course fixed by gauge 
invar iance. 

We need an ansatz for For the coordinate dependence, our ansatz will be similar to 
the one in refs. 10, 



u 



(P-l)/2 



(4.17) 



where /U runs over field theory directions, the ip{u)a are three spinor functions for which we 
must solve, and we extract a factor of u^^"^^/^ /~4 to make the Dirac equation look nice 
later (in the language of ref. ||2^ these factors will "remove the spin connection" from the 
equation of motion). 



Using our ansatz for the fermion in eq. ( 4.17 ) and the ansatz for the gauge field in eq. 
( |3.8D (from which we can recover eq. ( [3.12 ) simply by setting A], = 0), we find three Dirac 



^^Here the ± does not refer to the projectors n± of section^ but rather to the ± sign labehng the eigenvalues 
of the Dirac operator of S**^ in eq. (4.1C). 

^'^For a general (d-l-l)-dimensional AdS-Schwarzschild space, f{u) = 1 — u'^/u'^. 
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equations, 

= (^y^^du - ^ 7* + ^ka' - ^m) + ^ f ^^2, (4.18) 
= (^Vf^du - ^ 7* + - im) ^^2 - ^ 7* V'l + 7^ V-s, (4.19) 

= (v77"5« " ^ + " ^"') " ^^^""^ ^^-^"^ 

In what fohows we use the Lorentzian-signature from section ||, in which all the 
are Hermitian except for 7*, which will be anti-Hermitian:^'^ 

^-(t4)' ^-(t:). ^-(t:)' ^-(::)- 

(4.21) 

We will also use the projectors Hi 2, which in Lorentzian signature are defined as 

n« = i(l-(-l)"7"7V), (4.22) 

with Q = 1, 2. 

4.3.1 Normal Phase 

First consider the normal phase, where Al.{u) = 0. In that case, tps decouples from ipi and 
and its equation of motion becomes that of a free neutral fermion, as expected. We can 
then simplify the remaining two equations by taking linear combinations of them. Defining^^ 
^p± = tp2 i we find three decoupled equations, 

= (/77"5„ - ^ 7* + ^k^7' - ^m) V+ + ^ ^ 7* (4.23) 

= (//7"a, " ^ ^* + ~ l"^) ^- " ' ^^-^^^ 

= (^/77"5„ - ^ 7* + ikn' - lrr?j V'g (4.25) 

which are precisely the equations of motion for fermions (in the fundamental representation 
of the unbroken C/(l)3 C SU{2)) with charges q = =Fl,0. As in section ^, we emphasize 
that, because we consider a particular embedding of the Dirac equation into string theory, 
the allowed values of the mass and charge of the fermions are constrained by supersymmetry 
and T-duality. 



'We only need 7^* when P > 4. 



^^Here the ± index refers to Unear combinations of the worldvolume fermions that diagonahze the equations 
of motion when Al.{u) — 0. Nowhere in this subsection or the next do we use a ± index to refer to the 
projectors n± of section ^ or to the eigenvalues of the Dirac operator on 5''' of eq. (4.1C). 
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We will now follow appendix A of ref. [|11| to simplify the equation of motion further. 
First we rewrite the equation as^*^, 

(\/77"5n -^n?jiJ + iK^{u)ri^ = 0, (4.26) 



K^{u) = {-viu),ki), viu) = —{uj + qAtiu)), (4.27) 



1 

77 

where the index i runs over spatial directions, g = =Fl for 'ip± and q = for il^^. Notice that 
near the boundary, v{u) — t- w + g/i, so the frequency to is measured relative to {q times) the 
chemical potential. 

The system is rotationally invariant, so without loss of generality we can take only to 
be nonzero. (Obviously, this will not be the case in the superfluid phase, where rotational 
symmetry is broken.) The fermion's equation of motion then depends only on 7", 7* and 
7^, so the projectors Ha commute with the operator acting on the in eq. ( 4.26| ), hence 
the equations for (pa = I^ai^ decouple from one another. In terms of the (pa, the equation of 
motion becomes 

(^du + ^ma3^(l)a + ^ {-iv{u)a2 - {-If k^ai) (Pa = 0. (4.28) 

We thus obtain six decoupled equations, four for the (p±a and two for the (psa- 

Eq. (|4.2^ ) is almost identical to eq. (A14) of ref. The biggest difference is the 

function f{u), which for us is the f{u) of (4+l)-dimensional AdS-Schwarzschild and in ref. 
was the f{u) of (3+l)-dimensional AdS-Schwarzschild. Given that we will solve nearly 

identical equations of motion, we will obtain qualitatively similar finite-temperature results. 

As mentioned above, however, we cannot reach T = within the probe approximation, so 

we will not be able to reproduce the T = results of refs. ||9|, [l^, including in particular the 

influence of an emergent AdS2- 

4.3.2 Superfluid Phase 

In the solution corresponding to the superfluid phase, where Al.{u) is nonzero, we cannot 
write linear combinations of ^/Jl, ip2 and V'3 to diagonalize the system and produce three 
decoupled equations. To make comparison to the normal phase easier, we will again work 
with ijj± = ^2 i iipi, so that eqs. ( [1.18D -( p^ ) become 

= (y77"9„ - ^ 7* + ik^Y - V+ + i 7* + 4(^) 7^ ^3, (4.29) 

= (^^/fr^u - ^ 7* + ika' - im) V'- - i 7* V'- + 4(«) ^ i^3, (4.30) 

= (^^/fr^u - ^ 7* + ika' - im) Vs - ^ 7^ (V'+ + V'-) ■ (4.31) 



■"^Starting now, we will use the notation ip to refer to any of our three fermions, when we are making general 
statements. 
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Clearly the three fermions ^p± and ^3 couple to one another via a nonzero Al.{u). Here 
we have a concrete example of the couplings described in section (especially around eq. 



(4.32) 
(4.33) 
(4.34) 



( |2.32| )). We simplify the ip+ and ip- equations again by writing them as 
= (^y77"5„ - im^ ^+ + iK^{u)j^'i;+ + AI{u) 7" ^3, 

= (^y77"5„ - im^ V's + ii^,.(^)7^V'3 - I Aliu) 7" (V'+ + V- 



where is defined the same way as in eq. ( [4.27 ). Recall that tp^ has charge q = —1, ip- 
has charge q = +1 and ips has charge q = 0. 

Now we come to a big difference from the solution corresponding to the normal phase, 
at least for Dp-branes wrapping AdSp with P > 4. In the solution corresponding to the 
superfluid phase, rotational symmetry is broken from SO{P — 2) to SO{P — 3). Using the 
SO{P — 3) rotational symmetry, the most general momentum we can pick has nonzero kx 
and nonzero ky. The equations for xjj^ and ip- then depend on 7", 7*, 7^ and now also 7^, 
hence the H^ projectors no longer commute with the operators acting on the fermions in the 
equations of motion, and the equations for the (pa = Hq,^ will not decouple from each other. 
Upon acting with the projectors Hq,, the equations of motion become 



771 1 

du H 7= £73 + — {-iv{u)a2 - (-1)° k^ai) 



■'+a 



-^ai (-1) e 



(-l)"i<Ti(/<3a, (4.35a) 



TTl 1 

du H 7= 0-3 + ^ (-iw(n)(T2 - (-1)" k^ai) 

uVf Vf 



ky 



A. 



-^(-l)"2ai(/<3a, (4.35b) 



Tfi 1 

9« H 7= 0-3 + ^ (-i?j(n)cr2 - (-1)° k^ai) 



Ha 



1 



+ 



(4.35c) 



Here is antisymmetric with e^^ = +1. Notice that when ky is nonzero, the 0i and 02 
couple. 

Eq. ( 4.35| ) is the result for any Dp-brane extended along AdSp x . What will change 
from one Dp-brane to another are the allowed values of 771 and the solutions for A^(n) and 
A\{u). In the next section we specialize to the D5-brane extended along AdS/i x S"^ (P = 4 
and Q = 2), and to the massless worldvolume fermion. 
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5. Emergence of the p-wave Fermi surface 



5.1 Properties of the Spectral Function 

For the probe D5-brane worldvolume fermions of the last section, we solved the linearized 
equations of motion, eqs. ( [4.35 ), numerically, and used these solutions to extract the fermionic 



spectral functions. In this section we present a selection of our numerical results. 

We work with two D5-branes extended along (when T = 0) AdS^ x S"^ inside AdS^ x S^. 



As shown in section 4.1, we have many worldvolume fermions to choose from, with many 
different masses. In our numerical analysis we work exclusively with the massless worldvolume 
fermion. The dual operator is then the I = case of the mesino operator J^i written explicitly 



in eq. ( 4.13 ). These mesinos are valued in the adjoint of the SU{2) isospin symmetry, so 
we actually have three mesinos, , , and Jq, where the superscript denotes the charge 
under 1/(1)3. These are dual to the three fermions and in subsection 4.3.2 . 



In more detail, our procedure is as follows. We first choose the values of T and /x that 
we want, and solve for the background SU{2) gauge field functions A^{u) and Al.{u), as in 
section 3^. Plugging the gauge field solution into eqs. ( [4.35| ), we then solve for the bulk 



fermions. Near the horizon the fermions have the form of an in-going wave, eq. (g), with 
the a in that equation being a = for our system. When Al.{u) is nonzero, the fermions 
couple to one another, hence we employ the technique of section § to compute the retarded 
Green's function, which then gives us the spectral function, as we explain below. 

The normalization of our Green's functions is fixed by the normalization of the fermionic 



part of the D5-brane action, eq. (4.1). The normalization includes various numerical factors, 
and in particular depends on the normalization of the spinor x defined in eq. (|4.6D . 
We will omit the details, but we will mention that the normalization includes a factor of 
NfTD5 oc ^/XNfNc- In what follows we will rescale our Green's functions by the overall 
normalization. In other words, we will divide the action by the normalization factor, so that 
we obtain an effective AdS4^ Dirac action with a Lagrangian of the form i^p^. 

As explained in section |^, with three bulk fermions, the field theory retarded Green's 
function will be a 6 x 6 matrix. In the normal phase where the three bulk fermions decouple, 
the Green's function will be diagonal in isospin indices and in the subspaces defined by the 
III, 2 projectors defined in eq. (|2.12|) (see also eq. (|2.13|) ). Explicitly, the retarded Green's 
function will have the form 

[G%{u;, A:,, ky)] = diag {G^^, G^,, G% G^i, G^2, G^i) , (5-1) 

with A,B = 1, ... ,6, so that A = 1 corresponds to the components of the J-q mesino in 
the 112 subspace, A = 2 corresponds to J^q in the Hi subspace, ^ = 3 corresponds to the 
components of the J^q mesino in the 112 subspace, and so on. In the superfluid phase, due to 
the bulk couplings, for generic momenta all the off-diagonal elements become nonzero. 

The spectral function TZab is defined as the anti-Hermitian part of the retarded Green's 
function, 

TZab{uj, kx, ky) = i (^G%{uj, kx,ky) - G^^q{uj, kx, kyj^ . (5.2) 
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We define the spectral measure Tl{uj,kx,ky) as the trace over TZAsi^, kx, ky) (a trace over 
both flavor and spinor indices), 

TZ{u;, kx,ky) = trTZABi^, kx,ky). (5.3) 

Stability requires the eigenvalues of TZab, and hence both the diagonal elements of TZab 
and the spectral measure (the sum of the eigenvalues), to be strictly non-negative. Otherwise, 
if we perturb the medium with one of the operators appearing in the spectral function, the 
resulting excitation would experience negative energy dissipation into the medium, i.e. the 
excitation would extract energy from the medium, signaling an instability. The eigenvalues 
of the spectral function are a direct measure of the states of the theory that have an overlap 
with the relevant operators. The off-diagonal elements of TZab, however, need not obey any 
positivity requirement. Indeed, the off-diagonal elements are more similar to interference 
effects than to a measure of a density of states. We have confirmed numerically that all of 
our spectral functions obey the correct positivity requirements. 

We have confirmed that our numerical result for the spectral function obeys the following 
symmetries in both the normal and superfluid phases (in each case, any argument not shown 
is invariant): 

7^ll(-^J) = TZuiuj) , 7^22(-w) = 7^33(w), 7^55(-a;) = TZmi^) , (5.4) 
Tiiii-kx) ='R'22{kx), T^33{-kx) =T^u{kx), 'R'55{-kx) = T^eeikx) , (5.5) 
TZii{-uj,-ky) = TZu{uj,ky), 7e.22(-w, -A;^) = 7^33(a;, A;^), TZ55{-u;,-ky) = TZQQ{uj,ky) . (5.6) 

An example for frequency and momentum symmetries in the off-diagonal elements is 

nAB{-ky) = {-i)^+''nBA*{ky) , (5.7) 

which is also true for all diagonal elements, stating their invariance under ky — t- —ky. We 
have also confirmed that our numerical results obey 

TZab = TZba* , (5.8) 



which follows directly from the definition of the spectral function in eq. (^ 

Lastly, notice that because we study a fermionic operator of dimension 3/2, the retarded 
Green's function, and hence the spectral function and measure, are dimensionless. 

5.2 Numerical Results 

First we compute the spectral function for temperatures below Tc but in the normal (non- 
superfluid) phase, which we know is thermodynamically (disfavored. We do so for two reasons: 
first, to compare later to the superfiuid phase, and second, to reproduce some of the finite- 
temperature results of ref. as a check of our methods. In practical terms, we use the 
solution for the gauge field with A^{u) (from eq. ( 3.12| ) with P = 4) and zero A],(n). Figure § 
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shows two diagonal spectral function components, TZAAi^i kx, ky) with ^ = 1,5, as functions 
of kx/TrT, with Lv = ky = 0, for T < in the normal phase. 

Figure ^ (a) shows the component TZn, which corresponds to the components of the 
mesino J-^ in the 112 subspace, at temperature T = 0.61Tc. For the moment, our main 
point is that figure || (a) is qualitatively similar to figure 4 of ref. As we lower the 

temperature, the peaks in the figure move to larger momenta and additional peaks appear 
near zero momentum. In fact, in figure |2| (a), TZu along negative momenta (dashed blue 
curve) already displays a small bump near /c^/ttT = 0, which grows into a peak as we cool 
the system. Similar effects were observed in ref. Q, and were interpreted as the emergence 
of multiple Fermi surfaces at different momenta. Additionally, the spectral functions for our 
other charged fermions in the low-temperature normal phase are similar to those in refs. 
[P, (so we will not present them). 

Figure ^ (b) shows the component 7^55, which corresponds to the neutral operator Tq, 
in the 112 subspace, at temperatures T = Tc, 0.75Tc, and 0.61Tc. 7^55 is featureless here, but 
will not be so in the p-wave phase. Notice that 7^55 also does not change as the temperature 
decreases, or equivalently as the chemical potential increases, since when q = the chemical 
potential does not enter the relevant bulk fermion's equation of motion: see eq. ( [4.27 ). 




012345 012345 



kx/irT kx/irT 
(a) (b) 

Figure 2: Two spectral function components Ti-AAi'jJ, kx, ky), for A = 1, 5, plotted versus the rescaled 
momentum kx/i^T, with a; = fcj, = 0, at temperatures T <Tc but in the normal (non-superfluid) phase, 
i.e. in the thermodynamically disfavored phase, (a) TZn, corresponding to the fermionic operator 
with charge q = —1, at T = 0.61 Tc. The two curves are for positive momentum (solid red curve) 
and negative momentum —k^ (dashed blue curve), the latter case being equal to 7?.22(w, fc^,, /cy) at 
positive momentum due to the symmetries of the spectral function. Multiple peaks are visible, just 
as in ref. (b) TZ^^, corresponding to a fermionic operator with charge g = 0, at temperatures 

T = Tc (solid curve), T = 0.75 Tc, (grey dotted curve), T — 0.61 Tc (dashed curve). The curves remain 
coincident since changing the chemical potential does not affect the uncharged operator. 

Next we plot the essentially the same thing as in figure I, but now in the thermodynam- 
ically favored superfluid phase. More precisely, figure |3| shows two diagonal spectral function 
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components, TlAAi^,kx,ky) with j4 = 1,5, as a function of the rescaled momentum kx/irT, 
with uj = ky = 0, for T < in the superfluid phase, i.e. now with nonzero Al.{u). Figure 
^ (a) shows the same component of the spectral function as in figure |2| (a). Tin, again with 
T = 0.61Tc. Figure |^ (b) shows the same component of the spectral function as in figure ^ 
(b), 7^55, at the same temperatures T = Tc, 0.75Tc, and O.GlTc. 

The operator mixing is obvious in figure ^: the spectral function for a neutral fermion, 
7^55, develops a bump that grows into a small peak located at the same momentum as the 
peak in TZu, kx/nT = 3.87. In bulk terms, the coupling between (p^ and </>± in eq. ( [4.35 ) 
is allowing the peak in the charged fermions' spectral functions to "leak" into the spectral 
functions of the neutral fermions. That coupling is proportional to Al.(u), hence the peak 
should grow as the temperature decreases and Al.{u) grows, which is indeed what we see in 
figure ^ (b). The method we developed in section ^ for computing retarded Green's functions 
for coupled bulk fermions seems to work very well. 
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Figure 3: Two spectral function components ^ZAA{^^ k^, ky) plotted versus the rescaled momentum 
kx/irT with a; = fcy = in the superfluid phase, (a) Exactly the same spectral function as in figure 
H (a), TZii, at the same temperature T — 0.61 Tc, but now in the superfluid phase. The two curves 
correspond to positive momentum k^ (solid red curve) or negative momentum —k^ (dashed blue curve), 
(b) Exactly the same spectral function component as in ^(b), TZ^^, at the same temperatures Tc (solid 
curve), T — 0.75 Tc (grey dotted curve), T — 0.61 Tc (dashed curve). Here we see operator mixing: a 
feature develops in the neutral fermion's spectral function in the p-wave phase. A bump grows into a 
peak at the same momentum kx/i^T 3.87 as the peak in TZu. 



As another comparison between the normal and broken phases when T < Tc, we focus on 
the pole in the retarded Green's function Gii{uj,kx,ky), for the mesino with charge q = — 1 
in the 112 subspace, that is closest to the origin of the complex frequency plane, a; = 0, and 
follow the movement of the pole in the frequency plane as we change the momentum. 



''^Poles in the retarded Green's function are holographically equivalent to the bulk fermion's quasi-normal 



modes 



As explained in section ^ we can detect these quasi-normal modes from the vanishing of 



detP (e), where the matrix P is defined in eqs. (2.35) and (2.36). 
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Figure 4: The movement of the pole in the retarded Green's function Gii{uj,kx,ky), for a mesino 
of charge q = —I, that is closest to cj = 0, as a function of momentum, for T < Tc, in both the 
normal (thermodynamically disfavored) phase and superfluid (thermodynamically favored) phase, (a.) 
The movement of the pole in the normal phase at T = O.OlTc, for kx/i^T £ [2.58,12.58]. As the 
momentum increases, the pole moves from the lower right, near Rew/TrT = 0, toward the upper 
left. (The same applies for the following three figures.) The pole asymptotically approaches the real 
frequency axis, Imw/TrT = 0, as increases, (b) The movement of the pole in the normal phase at 
T — 0.48Tc, for kx/irT € [6.15, 16.15]. Here we see that the distance to the real axis does not decrease 
monotonically, but rather a local minimum develops at (Re w/ttT, Imw/TrT) = (—1.79, —7.21 x 10^^) 
when kx/irT = 8.15. (c.) The movement of the pole in the superfluid phase at T = 0.91Tc, for 
kx/T^T e [2.49, 12.29] and ky = 0. The movement is qualitatively similar to (a), (d) The movement 
of the pole in the superfluid phase at T = 0.48Tc, for kx/irT G [4.83, 8.98] and ky — 0. The movement 
is again qualitatively similar to (a), in particular, the distance to the real frequency axis does not 
develop a local minimum, in contrast to the normal phase result in (b). We see qualitatively similar 
behavior when we set kx — and increase ky. 



Figure ^ (a) shows the movement of the pole in the normal phase when T = 0.91Tc 
for values of hj.j'nT G [2.579, 12.580], which in the figure corresponds to starting at the point 
nearest Re[a;/7rT] = and moving toward the upper left. The pole asymptotically approaches 
the real frequency axis Imu/vrT = 0, as kx increases. 
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At a temperature T* ~ 0.6Tc, however, the distance to the real axis develops a local 
minimum. Figure ^ (b) shows the movement of the same pole as figure ^ (a) at a temperature 
T = 0.48Tc < r*, still in the normal phase. Here we see that the distance to the real 
frequency axis has a local minimum at (Re w/vrT, Im w/vrT) = (—1.79,-7.21 x 10^^) when 
kx/'KT = 8.15. Such behavior persists to lower temperatures, and indeed, the distance to the 
real frequency axis decreases. The lowest temperature we studied was T = 0.19 T^, where the 
local minimum occurred at (Re w/ttT, Im w/vrT) = (5.32,2.64 x IQ-^*^) when k^/TrT = 23.22. 

We seem to be seeing the emergence of a Fermi surface, which, as in ref. Q, would occur 
at T = when the pole would reach the origin of the complex frequency plane at some finite 
momentum kp, the Fermi momentum. Let us consider low temperature, and define k* as 
the value of momentum where the closest approach to the real frequency axis occurs. When 
T = 0, k* would be the Fermi momentum, k* = kp. At our lowest temperature, T = 0.19 Tc, 
k*/-KT = 23.22, and the closest approach to the real frequency axis occurs at an w* given by 
(Rea;*,Imw*) = (5.32, 2.64x 10"^^)7rr. Letting /c* play the role of /c^, we see behavior similar 
to the results of ref. Q: for small but nonzero temperature, as we change the momentum the 
frequency of the pole behaves as 

u-uj* ^ {k-k*y , (5.9) 

where our numerical results indicate that the exponent z = 1.00 it 0.01, and the spectral 
function behaves as 

nii^{k-k*y, (5.10) 

where our numerical results indicate that the exponent a = 2.0 it 0.1. In fact, these results 
are rather robust: we find the same z and a for many values of T < T*, and for several other 
poles. These results suggest that the low-temperature normal phase may not be a Landau 
Fermi liquid, which would have z = a = 1. 

Figures ^ (c) and (d) show the movement of the same pole, at the same temperatures, 
but in the superfluid phase. Figure ^ (c) shows the movement of the pole at T = 0.91Tc for 
kx/TrT G [2.49, 12.29] and ky = 0. Figure | (d) shows the movement of the pole at T = 0.48rc 
for kx/irT € [4.83, 8.98] and ky = 0. Unlike the normal phase result, here the distance to the 
real frequency axis does not develop a local minimum. In other words, here we do not see a 
Fermi surface emerge in the same fashion as in the normal phase. 

To see the emergence of the p-wave superfiuid Fermi surface, we study the spectral 
measure (largely following ref. [^8|), which as mentioned above, provides a direct measure of 
the density of states that have overlap with our fermionic operators. 

Our main results concern the evolution of the spectral measure TZ{u},kx,ky) as we cool 
the system through the superfluid phase transition. In the spectral measure we will see the 
breaking of rotational symmetry as we take T below Tc, and we will see the emergence of the 
Fermi surface as we approach T = 0, although in the probe limit we will not reach T = 0. 

^■^Our independent caluclations in the normal and superfluid phases yield the same position for the pole at 
T = Tcto within 0.1%. 
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Our results agree qualitatively with the T = results of ref . |18[] , in which the Fermi surface 
in the p-wave phase consists of isolated points. 

Figure ^ provides a road map for the evolution of the spectral measure as we lower the 
temperature. Here we set uj = 0, so we are studying fluctuations with zero energy above the 
chemical potential. Figure |5| (a) indicates the locations of peaks in the spectral measure, in 
the (fcj^/vrT, fcy/vrT) plane, with solid lines and the locations of small bumps as the dashed 
grey line. At T = Tc we see rotational symmetry: the black solid line indicates peaks for 
any values of (/caj/vrT, fcj^/vrT) on the black circle. At T = O.TTc the rotational symmetry is 
mildly broken: the green line is not a perfect circle. At T = 0.43Tc, sharp peaks only appear 
at isolated points on the axes, denoted by the red and blue dots (and also at the blue dot at 
the origin), while the dashed line indicates a small bump, rather than a sharp peak. Here we 
are clearly seeing the emergence of the Fermi surface at isolated points. 

To illustrate what the peaks and bumps look like, we choose a representative slice of the 
{kx / ttT , ky / ttT) plane, namely the line given by the polar angle (p = vr/S, which is drawn in 
figure [s] (a), and plot the (a; = 0) spectral measure versus the magnitude of the momentum 
1^1 = y''^^^~~'~~^ divided by vrT. In figure |5| (b) a distinct peak is visible for both T = Tc (black 
curve) and T = 0.7Tc (green curve), while for this generic (off-axis) value of (p the only feature 
of the spectral measure at T = 0.43Tc (dashed grey curve) is a small bump. On the axes 
{(p = 0,-7r/2), the picture is similar, except the bump becomes a sharp peak, corresponding to 
the red or blue dots in figure ^ (a) . 

To illustrate the evolution of the spectral measure in more detail, we present three- 
dimensional plots of TZ{u), kx, ky), for a; = 0, over the (kx/T^T, ky/nT) plane, for temperatures 
T = Tc (figure |(a)), T = 0.91Tc (figure | (b)), T = 0.69Tc (figure | (c)) and T = OATc 
(figure P (d)). In figure ^ (a) we see the peaks corresponding to the black circle in figure ^ 
(a). Clearly here the spectral measure is rotationally symmetric. When we cool the system 
to T = 0.69Tc (figure ^ (c)), we clearly see the emergence of five distinct peaks, two on the 
kx/irT axis, two on the ky/irT axis, and one at the origin. The circle of peaks corresponds 
to the green circle in figure ^ (a). When we further cool the system to T = OATc, the five 
peaks are still present, although the resolution of our three-dimensional plot is insufficient to 
resolve the two on the ky/nT axis away from the origin. 

Although these peaks have a much smaller footprint in the {kx/TrT, ky/irT) plane than 
the peaks on the kx/irT axis, they are much taller. The spectral measure TZ is of order 5 x 10^ 
at the peaks on the ky/irT axis, but only order 10^ at the peaks on the kx/irT axis, and order 
5 X 10^ at the peak at the origin. Apparently a large number of states are piling up at two 
precise locations on the ky/nT axis. 

For a bulk theory with SU{2) gauge fields and fermions in the fundamental representation, 
a combination of analytic and numerical results reveal that the a; = spectral measure at 
T = consists of two isolated points on the kx axis, located symmetrically about the origin 
[||]. For the same bulk theory but for fermions in the adjoint representation, as we consider 



here, the analytic arguments of ref. [18] indicate that the oj = spectral measure at T = 
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(a) 




Figure 5: The evolution of the spectral measure TZ{ijj,kx,ky) = trTZAsi'^, kx, ky) in the 
{kx / ttT , ky / ttT) plane at w = 0, as we lower the temperature, (a.) The position of peaks in the 
spectral measure TZ are indicated by the curves as we lower the temperature from Tc (black curve) 
through 0.7Tc (green curve) to 0.43rc (dashed grey curve). The T — 0.43Tc case exhibits a small 
bump rather than a sharp peak, except for points on the k^ and ky axes. We have indicated the 
bump with the grey dashed curve and the peaks with red and blue dots, including the blue dot at the 
origin, (b.) The spectral measure TZ plotted for a representative slice of the {kx/nT, ky/irT) plane 
(still with u — 0), namely along the line given by the polar angle (\) = 7r/8 drawn in (a.). We plot 
TZ versus the magnitude of the momentum |fc| = \ k'^+ ky divided by ttT, at T = Tc (black curve). 



T = 0.7Tc (green curve), and T = 0.43rc (dashed grey curve). 



should consist of three isolated points, one at the origin and two on the kx axis, located 
symmetrically about the origin. 



Therefore, what we see appears to be consistent with the results of ref. |18| for the 
structure of the spectral measure in the superfluid phase. The main obstacle to a direct 
comparison is the probe limit, which restricts us to finite temperatures: we do not know 
which peaks in our spectral measure persist to T = 0. Nevertheless, given that we see the 
three peaks on the k^/nT that we generically expect, and that the two peaks on the ky/nT 
axis have a shrinking footprint as we cool the system, we have good reason to believe that 
the results of ref. may apply to our system. To answer the question fully requires 

computing the back-reaction of the D5-branes in the bulk. Whether that produces a domain- 
wall geometry of the kind used in ref. [^8| is not guaranteed. 

We begin to explore the u dependence of the spectral measure Tl{uj,kx,ky) in figure ^. 
Figure |7| (a) simply reproduces the T = 0.4Tc part of figure |5| (a), where ui = 0. Figure (b) 
shows the spectral function along three lines in the {kx/TrT, ky/nT) plane: along the positive 
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(a) T = r„ = 



(b) T = 0.91 Tc, Lo^O 








(c) r = 0.69 Tc, oj = 



(d) T = 0.43 Tc, a; = 



Figure 6: Three-dimensional plots of the spectral measure Tl{uj, k^, ky) in the superfluid phase over 
the {kx/'!TT,ky/nT) plane at zero frequency, lu — 0, and for distinct temperatures T < Tc- (a) The 
T = Tc case, which is clearly rotationally invariant. The peaks correspond to the black circle in figure |^ 
(a), (b) The T — 0.91Tc case, where the spectral measure does not yet display any dramatic breaking 
of rotational symmetry. Notice the peak a the origin of the {k^/TrT, ky/nT) plane, (c) The T = 0.69Tc 
case, where the breaking of rotational symmetry is obvious. We see that the "cylinder" of (a) breaks 
into five distinct peaks on the k^/nT and ky/irT axes. The peaks correspond to the green circle in 



figure 5 
figure 5 



(a), (d) The T — QATc case, which still has five peaks, labeled by the red and blue dots in 
(a). The peaks along the ky/irT axis are too narrow to appear in the three-dimensional plot 



with the resolution we use. 



kx/T^T axis (red curve), along the positive ky/nT axis (blue curve), and along the line given 
by the polar angle (/> = vr/S in figure I (a). Here we see explicitly the difference in widths of 
the peaks on the k^/irT and ky/nT axes (the red and blue peaks). 

To explore the uj dependence, we choose a few representative points in the {kx/nT, ky/nT) 
plane and, for each point, plot the spectral measure versus uj. For these points, we consider 
not only T = OATc, as in figure ^ (a), but also the slightly higher temperature T = 0.55Tc, 
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Figure 7: (a) Peaks in the w = spectral measure in the plane at T = 0.4Tc. The labeling 

is the same as figure ^ (a), (b) The w = spectral measure along three lines in (a): along the |y axis 
(red curve) and the ^ axis (blue curve), so the peaks at nonzero momenta correspond to (3) and (4) 
in (a), respectively, and along the polar angle 4> = ""/S (dashed grey curve), drawn in figure ^ (a), so 
the bump corresponds to where </) = 7r/8 intersects the dashed grey line in (a). In each case the peak 
at zero momentum corresponds to (0) in (a). (1) 7?. as a function of lo/ttT for the point labeled (1) in 
(a). The dotted and solid lines arc for T — O.SSTc and T = QATc, respectively. (2) and (3) show the 
same thing for the corresponding points (2) and (3) in (a). 
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in order to study the behavior as we cool the system. 

Our points are similar to those in figure 8 of ref. [18|, where the same quantities were 
plotted (for the slightly different system of ref. the spectral measure versus uj for fixed 

kx and ky. We will thus compare our results to those of figure 8 of ref. [^] along the way. 

In figure |7| (1) we plot the ui dependence of 7^(a;, kx, ky) for the point (1) labeled in figure 
^ (a). The dotted line is for T = O.SSTc and the solid line is for T = 0.40rc. Here we see that 
as we cool the system, a small gap (a depletion of states) opens near w = 0, while a sharp 
peak emerges near cj/vrT 1.7. Such behavior at least appears to be approaching that of 
figure 8 (1) in ref. [|l8|, where a genuine gap (zero states) appeared at w = 0. 

In figure |7| (2) we plot the uj dependence of lZ{uj, k^, ky) for the point (2) labeled in figure 
(a). The dotted line is for T = 0.55Tc and the solid line is for T = 0.40Tc. Here we see 
that a sharp peak near oj/ttT w 1.5 when T = 0.55Tc shrinks and begins moving toward 
a; = as we lower the temperature to T = OATc. Moreover, the small peak near lv/ttT ~ 3.5 
when T = 0.55rc grows much sharper at T = OATc- As in ref. ||T^, here we seem to see the 
emergence of the well-known "peak-dip-hump" shape, with the peak being at uj/ttT « 3.5, 
the dip being at oj/ttT ^ 5.2, and the hump being at oj/ttT 6.5. 

In figure 8 (2) of ref. [18|, a gap was present in the spectral measure for small frequencies, 
except for a single genuine delta-function peak at finite frequency, and at larger frequency 
a continuum of states appears (the "hump"). As argued in ref. [|l^], at finite temperature 
the delta-function peak will acquire a finite width and merge with the hump, producing the 
peak-dip-hump. Our spectral measure appears to be approaching the form of the spectral 
measure in figure 8 (2) of ref. [jl8| (with the usual caveat that we cannot actually reach 

r = o). 

In figure |7| (3) we plot the u dependence of TZ{uj, kx, ky) for the point (3) labeled in figure 
(a), which is sitting right on top of the peak on the positive kx/T^T axis. The dotted line is 
for T = 0.55rc and the solid line is for T = 0.40Tc. At the higher temperature (the dotted 
line), the primary feature is the peak near uj/ttT ~ 1.2, which moves toward a; = and also 
shrinks (the peak is lower) as we lower the temperature, becoming the peak in the solid line. 
Assuming such a trend continues, our results would be consistent with figure 8 (3) of ref. |18], 
where, sitting right on top of the peak on the kx axis, the spectral measure went to a finite 
constant at a; = 0. 

Figure ^ shows TZ{u}, kx, ky) versus uj/ttT for the points (0) and (4) in figure (a), sitting 
right on top of the peaks at the origin and on the positive ky/nT axis, respectively. Here we 
use only T = OATc. We clearly see a gap developing at low frequency in both cases. Such 
behavior is similar to the gap that develops in the spectral function of vector fiuctuations 

|3^, which is immediately related (via a Kubo formula) to a gap in the conductivity. 
Whether these two gaps are related is unclear, but deserves further study. 

Finally, to explore further the u dependence of 7^(a;, kx, ky), we do not restrict to points in 
the {kx/TrT, ky/irT) plane, but rather restrict to a single nonzero value of frequency, uj/ttT = 
0.25, and plot the spectral me asure over the entire {kx/irT, ky/nT) plane. The result appears 
in figure where (a) is for T = Tc, (b) is for T = 0.91Tc, (c) is for T = 0.54Tc and (d) is 
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Figure 8: The spectral measure TZ{uj,kx,ky) as a function of uj/ttT for {kx/TrTjky/irT) values 
corresponding to points (0) (dotted blue curve) and (4) (solid blue curve) in figure ^ (a). These points 
sit right on top of the peaks at the origin and on the positive ky/nT axis, respectively. Here T = OATc- 



for T = 0.4Tc. We see a number of differences from the w = case of figure ^ At the 
transition, T = Tc, the spectral measure is again rotationally symmetric, but now with two 
concentric "cylinders." By the time we cool the system to T = 0.54Tc, we see a number of 
peaks clumped near the momentum axes. If we compare the spectral measure at T = 0.4Tc 
at a; = and co/irT = 0.25, figures |^ (d) and |9| (d), respectively, then we see that the a; = 
peaks on the ky/nT axis each split into a number of peaks at oj/ttT = 0.25 which then move 
apart along the kx/i^T axis. 

To summarize, in the superfluid phase we observe isolated peaks in the spectral measure, 
whose locations appear to be consistent with the p-wave nature of the condensate, as well as 
with the fact that our fermions are in the adjoint representation of SU{2). Such a structure 
suggests nodes in the energy gap on the normal-phase Fermi surface. We plan to investigate 
this further in the near future, for instance by identifying the appropriate Dirac cones and by 
studying the spectral measure for fixed values of as well as for further values of fixed kx, 



6. Conclusions 

The main results of this paper are: 
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(c) T = 0.54 Tc, oj/ttT = 0.25 (d) T = 0.43 T^, oj/ttT = 0.25 

Figure 9: Three-dimensional plots of the spectral measure 7?,(lli, kx,ky) in the superfluid phase over 
the {kx/irT, ky/nT) plane at lo/ttT ~ 0.25, and for distinct temperatures T < Tc. The four plots are 
arranged in a similar fashion as in figure |[ (a) is the T = Tc case, (b) is the T — O.OlTc case, (c) 
is the T = 0.54Tc case, and (d) is the T = OATc case. The main differences from figure ^ are that 
the "cylinder" we saw in figure ^ (a) is now two concentric cylinders, and the peaks along the ky/nT 
axis each split into multiple peaks which then move away from one another along the k^/T^T axis as 
we cool the system. 



• We described an efficient algorithm for computing retarded Green's functions from bulk 
fermions that couple to one another, 

• We performed holographic renormalization for fermions in AdS, 

• We determined a simple embedding into string theory of a charged fermion in AdS, 
using probe D-branes, 

• Using the above ingredients, we numericahy computed fermionic spectral functions in 
p-wave superfluid states, with emphasis on the emergence of the isolated points of the 
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Fermi surface as we cooled the system through the p-wave superfluid transition. 

Using the machinery that we have developed, we can think of a number of directions for 
future research. 

Introducing a magnetic field on the worldvolume of a probe Dp-brane is straightforward 
to do, and has many interesting consequences |34, 43, ^ 48, ^ 90, |9^, 94, 95]. 



Our embedding of p-wave superfluids, and of charged fermions, into string theory may be 
a useful arena in which to study the effects of magnetic fields both on holographic p-wave 
superfluids, along the lines of ref. |34|, and on holographic Fermi surfaces, along the lines of 
refs. m, m, ||. 

Typically the mesinos carry some R-charge, so another avenue to study holographic Fermi 
surfaces would be to embed probe Dp-branes into background geometries corresponding to 
field theory states with finite R-charge density. In those cases the T = limit is accessible 
within the probe approximation. These geometries will generically be charged dilatonic black 
hole solutions. The most attractive of such solutions, for condensed matter applications, may 
be the solution discussed in ref. |100|] , which produced in the field theory specific heats linear 
in temperature. 

We focused on probe D5-branes extended along AdS^ x 5^ inside AdS^ x S^, and worked 
only with the massless worldvolume fermion. Other Dp-branes, and worldvolume fermions of 
other masses, deserve further study. 

Lastly, knowing the T = 0, finite isospin density ground state of the D3/D5 theory would 
be extremely useful. In bulk terms, the problem is to find a solution of type IIB supergravity 
representing fully back-reacted D5-branes with non-trivial worldvolume gauge fields. 

We plan to investigate these and related issues in the future. 
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A. Holographic Renormalization of Fermions in AdS 

In this appendix we perform the holographic renormalization of a Dirac fermion in spaces 
that asymptotically approach AdSd+i- We will work in Euclidean-signature AdSd+i, unless 
stated otherwise. We assume the bulk metric asymptotically approaches the metric of eq. 

dtp 1 

ds^ = QAB dx^dx^ = — TT H — dx^dx^ , (A.l) 
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but, unlike section ^ we leave the number of spatial dimensions d unconstrained. 

Our results will only be valid for spaces whose metrics asymptotically approach the 
AdSfi-\-i metric. Various generalizations are possible. One generalization would involve spaces 
for which surfaces of fixed u are not simply flat space, that is, spaces for which 6ij — t- gij{x) in 
eq. ( [A.ip , where gij{x) has non-trivial curvature. Another generalization would be to study 
spaces that are not asymptotically AdSd+i, for example, spaces that asymptotically approach 
the near-horizon geometry of Dp-branes, with p < 5 but p ^ 3. We also leave an analysis of 
the Lorentzian-signature case, along the lines of ref. |101| , |102| , |103| ], for the future. 

As mentioned in section |2|, the form of the counterterms will be fixed by symmetries (in 
particular Lorentz invariance on the u = e surface) and by the requirement that they cancel 
divergences of the on-shell action. Moreover, for fermions we have an additional constraint. 



The counterterms can be built only from (defined in eq. 2.13 ), because that is held fixed 
under variations, and hence does not spoil the stationarity of the action. 

Otherwise, the procedure for fermions very closely parallels that for scalars, which was 
reviewed for example in ref. [61]. Our main results are the counterterms in eqs. ( |A37|) and 
( |A.34 ) and the renormalized on-shell actions in eqs. ( lOsI) , ([091) , and (lAlCp . 



A.l Solving the Equation of Motion 

We begin with the Dirac action (plus boundary terms). 



bdyi 



(A.2) 



We will work with a single Fourier mode, so we let ^' — )■ e^^^ where, without loss of 
generality, we have chosen the momentum to point in the x direction. The equation of 
motion is then eq. (|2.8|) , 



^' = 0. 



(A.3) 



As in section |2| we define = ^ (1 it 7") ^ so that 7"^^ = ±^±. We will refer to and 
as having opposite chirality (even when d + 1 is odd). With we obtain the same 
decoupled second-order equations as in eq. (|2.18|) , 

= 0. (A.4) 
We can solve eq. ([A.4|) exactly. The form of the solution depends on the value of m. 



u 



2 d 
-m lb m H — - H — 



First, suppose m is not a half-integer. The solution of eq. (A.4) is then 



(A.5) 
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where the J functions are Bessel functions of the first kind^^ and C^{k) and C^{k) are spinors 
of the same chirahty^^ as which may depend on k, as indicated. 

When m is half-integer, m it ^ is an integer and hence the order of the Bessel functions 
is an integer. In that case the two Bessel functions are not linearly independent: if the order 
n of Jn{x) is an integer, then Jnix) obeys the special identity J_„(x) = (— l)'^J„(x), so in 
that case Jn{x) and J_„(2;) are obviously linearly related. When m is half-integer, we must 
introduce a Bessel function of the second kind, Yn{x), which is linearly independent from 
Jn{x). When m is half-integer, the solution thus becomes 



d+i 



Ct{k)J^^.{^n^u) + Ct{k)Y„^^.{^r^u) . (A.6) 



Notice that if the order n of Yn{x) is an integer, then Yn{x) obeys a relation similar to that 
for J„(x), namely Y_n{x) = {-l^Ynix). 

The key difference between the Bessel functions of the first and second kinds that will be 
important for us is that the small-x expansion of Jn{x) involves only powers of x, whereas 
the expansion of Yn{x) involves both powers of x and logarithms of x. Indeed, Jnix) has a 
series expression (for any n), 

from which we can immediately read the expansion for small x. Yn{x) has a series expression, 
when n is a non- negative integer, 

k=0 

-;(l)"g^Fwi«"+'=+>'+*<'+'>i(l) • 



'''^Notice also that since the arguments of the Bessel functions are imaginary, we could equally well write 
them using the modified Bessel functions In{x) and Knix). For our asymptotic analysis, the distinction is 
not important. If the space was pure AdSd+i (not just asymptotically AdSd+i), however, then eq. ([A.5|) is 
the solution for all u, and regularity in the bulk would force us to discard the In{x) solution. Indeed, we will 
discuss this later in the appendix. 

'^''in the rest of the paper, 1 and 2 subscripts referred to the projectors ni,2. We never use these projectors 
in this appendix. Here, the subscripts 1 and 2 on C^{k) and C^{k) simply distinguish the two linearly 
independent solutions of eq. (A. 4). 



^^The Wronskian 

sin nn 



M^(J„(a;),J_„(x)) = 



TVX 

clearly vanishes when n is integer, indicating linear dependence. The Wronskian 

2 



W{Mx),Y„{x)) 

nx 



for any n. 
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where tl'{x) = digamma function. We can obtain the series form of simply 

from yL„(x) = (— The two features to notice are that the logarithmic terms in 
Yn{x) are simply of the form Jn{x) logx, and Yn{x) also includes terms with inverse powers 
of X, from the x~"' multiplying the first sum. 

Using the series expressions for the Bessel functions, we can rewrite the solutions in a 
form useful for our purposes (and introduce some notation to keep expressions simple and 
compact). First, consider m not-half-integer. We write 



± 



d+l 
U 2 

d+l 
U 2 



^=F2 

- 1 



(A:)n"T2 {l + s^{u,k)) +4(k)u 



- mT 



[l + s^{u,k)) 



= cf{k) ni+'"^5+5 (1 + st(u, k)) + c^{k) (1 + k)) , (A.9) 

where in the second line we have absorbed various factors into C^{k) and C^{k), which we 
then relabeled as cf{k) and ^^{k), and we have defined the series 



st{u,k) = £a±(m) i-k^u^^, af{m) 



j!22j- r(j + l+ (m^i))' 



(A.IO) 



=F \). We have 



and s^{u,k) and b^{m) are defined similarly, but with (m =F ^) — {m 
isolated the leading powers of u near the boundary n — )• 0: the sums s^{u,k) and s^{u,k) 
involve sub-leading powers of u (each sum starts at order n^). Written separately (to facilitate 
counting powers of u), ^± are 



^+ = c+ik) (1 + s+(u, k)) + 4ik) u^-"'+^ (1 + s+(n, k)) 



(k) u 



1 + (n, A;)) + C2 (k) u'- 



-{u,k)) 



(A.ll) 



The coefficients in eq. ( [A.ll ) are actually not independent. If we return to eq. ( |A.3| ) and 
insert ^ = ^+ + using our solutions for then we find, upon generalizing kj^ — t- ^ 
and collecting powers^^ of u, 



[{-2m + 1) c^{k) + i^C2ik)] ui~ 



■m+l 



+ [-{2m + 1) c^{k) + i¥4{k)] n2+™+i + O [u2~'^+'^ 
We thus conclude that 



2m + 1 



(A.12) 



(A.13) 



^^When m is an integer, another term of order it 2 +'"+'^ appears in eq. (A.12), with a coefficient proportional 
to b^{m)c2 (k) + b^{m)i]f! C2 (k) . This term has opposite chirality from the term shown in eq. (A.12), hence we 
set them to zero independently. Using the definition of the b^{m), we then recover exactly the same relation 
between C2'(fc) and c^(fc) as in eq. (A.13), so we obtain no new information. 



- 54 - 



When m is half-integer, we can plug the series expressions for the Bessel functions in eqs. 
into the solutions for and rearranging various terms, we can write 



O) and (I 



d+l 

U 2 



d+l 
U 2 



cf(A;)n"^^ (Inn) (l + s±(m, A:)) + ^(fc) ^-('"^s) {l + s^{u,k)) 
= cf(A;)ui+'"T3+l (Inn) (l + s±(n,A;)) +4(A;)nt-"^±^+3 (l + s^(n. A:)) , 

where in the second line we have absorbed various factors into C^{k) and C^{k), which we 
then relabeled as c^(fe) and c^{k). Notice in particular that we exchanged the indices 1 and 2. 
Recall that we are using units in which the radius of AdS^+i is equal to one. The arguments 
of the logarithms include factors of the AdSd+i radius to render them dimensionless. The 
sums s^{u,k) are defined in eq. ( [A. 10 ), while the sums s^{u,k) are 



s^{u,k) = ^df{m,k){-k 



2Vn2i, 



(A.14) 



where the coefficients df{m,k) are not particularly illuminating to see, so we will not write 



them. They can be derived straightforwardly from eqs. ( |A.7| ) and ( A. 8 ). That derivation 
also shows that they can depend not only on m but also on k, as indicated. We have defined 
our notation to isolate the leading powers of u near the boundary n — >• 0, and also to show 
that the solutions are identical in form to those for the m- not-half- integer case in eq. ( A. 9)) , 
except for an extra logarithmic factor in the c^{k) terms. Written separately (to facilitate 
counting powers of u), are 

^+ = c+(A:)nt+"^ (Inn) (l + s+(n. A;)) + c+(A:) ni-"+i (l + s+(n, A;)) 



Ci (k) u 



(Inn) (l + s-(n,A;)) +C2(A:)n2-'"(l + s^(n,A;)) . (A.15) 



Here again, the coefficients in eq. ( [A.15 ) are not independent. Inserting the solutions for ^'-i- 
into the Dirac equation, eq. (|A.3D , and generalizing kj^ — t- J/, produces 



[-{2m + l)c^{k) + il/ct{k)] nf +™+Mnn + O (^nf 



+m+2 



Inn 



+ (-2m + 1) c^{k) U2 
+c+{k)ui+"' + 4{k)J 



l + s+{u,k)) +i^c^{k)u2-'^+^ ^i + s^{u,k)) 
[n9„s|(n. A;)] + O 



-m+2 



(A.16) 



where in the first line we have indicated the sub-leading order of terms involving powers and 
the logarithm of n, while in the third line we have indicated the sub-leading order of terms 



involving just powers of n. When m ^ 1/2, the vanishing of the u' 
terms in eq. (A. IE) requires 

1 1 



\+m+l 



Inn and n2 



-m+l 



2m + I 



2m 



(A.17) 
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which are identical to what we found in the m non-half-integer case, eq. (|Al^ ). In what 
follows, we will also need to write cf{k) in terms of c^(A:). We thus turn to higher-order terms 
in eq. ( A. 16 ). When m / 1/2, the u^^"^ term will always be the same power of u as ii2~™+^ 
times a particular term in the summations s'^{u,k), s^{u,k) and uduS^{u,k). Recall that 
these summations involve powers of v? . Some term in the summations, say the j'th term, will 
have j = m — i. Recalling the definition of the summations in eq. ( [A. 14 ), we find from the 
second and third lines in eq. (A. 16) 



4{k) 



4{k){-2m + l)d'l{m,k) + il/c2 ik)d~ {m,k) + 4{k)dUm,k) 2j 



-i^C2 {k)d ^^{m,k) {-k"^) 



(A.18) 



where we plugged j 



m- 



^ into the first line, so that the two c^(/c) terms canceled, producing 



the second line. The factor of 2j = 2m — 1 in the first line is crucial for the cancelation: it 
comes from the udu acting on sj'(m, k) in the third line of eq. (A. 16). The m = 1/2 story 
is similar: ii2~™+i and ii2+™ become the same power 112 + 2^ so that, from the second and 
third lines in eq. ( A.16| ), we immediately find 

ct{k) = -ili^c^ik). 



(A.19) 



In the following subsections we will restrict to positive values of m unless stated otherwise. 
We can recover results for negative m as follows. For m not-half- integer, if m < 0, then we 
can insert m = —\m\ in eq. (A. 5) to obtain 



d+l 
U 2 



u— [Ct{k) J_(|„|±i)(V^n) + Cf{k) J|^|±i (v^n) 



(A.20) 



To obtain results for negative m, we can work with positive m and then in all formulas take 
m — )• \m\ and exchange C^{k) o C2^(/c), which means 

cf{k)^4{k), 4{k)^cf{k), 



For m half-integer, if m < 0, then we can insert m 



m not-half- integer. 
\m\ in eq. ([A.6|) to obtain 



(A.21) 



d+l 
U 2 



Ct{k)J^^UV^u) + Ci{k)Y.{V^u) 



d+l 
U 2 



Cf{k) J_(|„|ii)(\/^^x) + Ct{k) ll(|^|±i)(v^n) 



Ct{k)J\^\^i{^r^u) + Ct{k)Y^^\^i_{^/^u) 

where in the third equality we have used J-n{x) = (— l)"Jn(a;) and Y^nix) = (— l)"l^(a;). 
To recover results for negative m, we can work with positive m and then in all formulas take 
m ^ |m|, Cf{k) (-l)l™l^icf (A;), and C^{k) {k) , which means 

cf (A;) ^ (-l)l"l^^cf (A;), ^(A;) ^ (-1)I™I^^4(A;), m half-integer. (A.22) 
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A. 2 Determining the Counterterms 

As reviewed in section tlie AdS/CFT dictionary equates the exponential of (minus) the on- 
shell supergravity action with the generating functional of field theory correlation functions. 
For the action in eq. ([A. 21) , clearly the bulk term vanishes when evaluated on a solution. The 
only nonzero contribution to the on-shell action comes from the boundary terms. We will 
split Sbdy into two terms, 

Sbdy = Svar + SqT, (A. 23) 

where Syar stice terms required for the variational principle to be well-posed |56, 33 1 and Sct 



includes the counterterms that will cancel any divergences. As shown in refs. [56, 63 



(A.24) 



where the integration is over only the u = e surface, ^77 is the determinant of the induced 
metric on the u = e surface, which for us is simply = e~'^, and are evaluated at m = e. 
We will plug the solutions for ^>± into S'^ar and isolate any terms that diverge as e — )• 0. We 
will then introduce into Sct any terms necessary to cancel the divergences and hence render 
the on-shell action finite. The terms in Sct must respect the symmetries of the on-shell Svar 
and must be built only from to preserve stationarity of the action, as explained above. 



For m not half- integer (and positive), if we plug the solutions for ^± from eq. ( A. 11 ) 
into eq. (|A.24| ), we find^^ 



Sit 



c{ C2 (1 + (e, k)) + c+c^ e 



(1 + (6, /C)) 



(A.25) 



d+2m+l 



(l + /a+a- (e,A:)) + c 



Co e 



d-2m+l 



2 ^2 



(1 + A+,- {e,k)) 



where we have introduced one more piece of notation: we have defined 

fa+a-{e,k) = st{e,k) + s~{e,k) + st{e,k)s~{e,k), 

and similarly for fh+b~{^->k)i fa+b-i^^k), and fa-b+i^^k), all of which are summations in 
powers of starting with e^. We now ask what happens to each of the terms in the brackets 
in eq. ( A.25| ) as e — )• 0. The first term clearly remains finite (the e'^'s cancel). The second 
and third terms vanish as and g^m+i^ respectively. The fate of the fourth term depends on 
m. If m < 1/2, then the fourth term vanishes as g-^m+i^ done: no divergences 

appear and no counterterms are necessary. If m > 1/2, however, then the fourth term may 
have one or more divergent terms, with the leading divergence going as g-2m+i^ must 
add counterterms. Starting now, we will assume m > 1/2. Using eq. ( A.13| ), we can write 
the on-shell Svar as 



Sv 



1 



2m - 1 



d-2m+l 



{1 + fb+b- {e,k)) + 0{e 



A+2 



. (A.26) 



''Starting now, we drop the k dependence in the coefficients: (fc) — >■ and cj(fe) —J- Cj 
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Having isolated the diverg ences in Syan now must write an S(jx that obeys the 
constraints mentioned in section ^: it must cancel all divergences of, while respecting all 
symmetries of, the on-shell Syar, and must be built only from the boundary value of The 
ScT that does the job is 

/oo ^ ^ oo 

d'^x^ aj(m)^_^,q^'_ = / (fx^ ^ei+2Jaj(m)'l'_^n%_, (A.27) 
i=o i=o 



where (^^ = egl (the power of e comes from the inverse vielbein evaluated at u = e) and 0^ 
is some power j of the scalar Laplacian on the u = e surface (so the derivatives act only 



in field theory directions), which in our case is simply = e d . When we take ^ — t- e* 
the counterterms become 



j=0 



(A.28) 



so that, plugging in the solutions for we find 



1 °° 

ScT = + d'^x^ j;6^+2ia,(m) 



d+2m+2 



d+l 



d+l 



d—2m 



^^il/{-kyc^{l + f,-,-ie,k)) 
^^^l/{-kyc^{l + h-,-{e,k)) 
c^tJ/i-kyc^il + h-,-ie,k)) 



(A.29) 



Of the terms in brackets, the term that is potentially divergent when e — t- goes as g'^-2m 

of the other terms vanish when e — t- 0. We fix the coefficients aj{m) by demanding that the 

divergent terms in eqs. ( A.26| ) and ( |A.29 ) cancel each other, which means that the quantity 



1 



2m - 1 



(l + A+6-(e,^)) + 



oo 

Ya,{m){-e'ky {1 + f,-,-{e,k)) 

j=0 



(A.30) 



must vanish order- by-order in — e^fc^, up to order e^™~^. We immediately see that ao{m) = 
— 2m-i ; and we can write a formal recursive solution for all the other aj{m) (here we drop 
the dependence on m for notational clarity, so aj{m) — )■ aj, etc.). 



a,- 



(2m - 1) 



Kj 



2 b: 



(A.31) 



k=l 
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where we define b^{m) = if j < 0. Tlie first four aj{ni) are 
1 K + K 

as = ^ - 02 - ai (26^ + &r ) " "o (263 + 26^6^) . 



Plugging in the explicit forms for the b^{m) (which are just the a^{m) from eq. ( A.10| ) with 
m =F ^ — — ("i T ^)), we can write these explicitly: 

1 

ao(m) : 



2m - 1' 

ai(m,) = 



(2m- l)2(2m-3)^ 

a2(m) 



(2m- l)3(2m-3)(2m-5) 
17 - 10m 

a3(m) 



(2m - l)4(2m - 3)2(2m - 5)(2m - 7) ' 



We have thus determined the counterterms when m is not half-integer. 

We now consider half- integer m, in which case the solutions for ^'-i- appear in eq. ( [A. 15 ) 



Oncer again, to determine the divergences, we plug the solutions for ^± into S^ar-, with the 
result 

Svar = j d'^X^i!+^^ 

= I d'^xl [c+C2 e'^ (Ine) [1 + f^^,.{e,k)) + c+c^ e^+^ (Ine) (1 + (e, A;)) 

where fa+d~{^,^), fa~d+{(^,k) and J^+rf- (e, A;) are defined simil arly to f^+a- {^,k). We now 
ask what happens to each of the terms in brackets in eq. ( A.32| ) when e — )■ 0. The first term 



diverges as Ine. The second and third terms vanish as Ine and e^m+i (ine)^, respectively. 
The fate of the fourth term depends on m. If m = 1/2, then the fourth term is finite. The 
first term still diverges, however, so in this case we need counterterms (in contrast to the 
m- not-half-integer cases). If m > 1/2, then the fourth term may have one or more power-law 
divergent terms, with the leading divergence going as g-2m+i Notice, however, that the 
fourth term will also always produce something finite, since the sum fd+d~ (^j ^) will always 
have a term that goes as e^m-i -^^j^ich will cancel the g-2m+i_ \y]2en m 7^ 1/2, we can use 



eqs. ( A. 17 ) and ( A.18 ) to rewrite the S^ar in eq. ( |A.32| ) as 
1 



cic-e'' (Ine) (1 + (e, A:)) 



di X ^ C2 iiJj^ C2 



^ e'^lne(l + W(e,A;)) (A.33) 



+- 



1 



d-2ni+l 



2m -1 



2\ rra- 



The m = 1/2 case is similar, except in the third hne i{m,k) {—k"^) 
fourth hne remains identical to the first line. 
We now include logarithmic terms in Sct, 

/oo 
d^x [aj (m, A;) + {m, k) In e) ^ 
j=0 

Plugging the solutions for into ScTi we find 

00 

Sct = + d'^x^ J^e^+^J {aj{m,k) + (3j{m,k) Ine) 

+6^+1 (Ine)i -c-il^{-k^yc^{l + f,-,-{e,k)) 
+6^+1 (Ine)i c2i^(-A;2)J-cr(l + /,-.-(e,fe)) 



1, and the 



(A.34) 



d—2m 



(In ef c^i]/{-ky c, {1 + h-d-{e,k)) 



d-2m+l 



/ '^'^^ ~d X] («j("^> ^) + k) In e) e^-' e' 



C2-z^(-fc2)^ C2 (1 + (e, fc)) + O (e2™+i Ine) 



(A.35) 



Of the terms in brackets in the first equality, only the fourth term is potentially divergent 
when e — 7- 0. All the other terms vanish when e — )• 0. In the second equality we have thus 
isolated the potentially divergent terms. We now compare eq. ( [A.33 ) and eq. ( [A.35 ), and 
adjust the coefficients aj{m,k) and /3j{m,k) such that all divergences cancel. 

We first consider the logarithmically-divergent terms. For these we must adjust the 
coefficients f3j{m,k) such that 



' i=o 



(A.36) 



We thus conclude that all (3j{m,k) are zero except for one, 1 (m, A;) = d^^_i{m, k). If 
m = 1/2, the result is that only /3o(m, A;) = 1 is nonzero. 
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Next we consider the power-law divergences, which appear when m ^ 1/2. For these we 
adjust the coefficients aj{m,k) such that the quantity 

_ oo 

^— (1 + fd+d- (e, k)) + <^M^ k){-e^ky (1 + (e, k)) (A.37) 

vanishes order-by-order in — e^/c^, up to order g2m-3^ guarantee that all divergences can- 
cel. The aj{m,k) are thus straightforward to obtain, indeed, they are identical in form to 



the aj{m,k) in eq. (A. 31), but with b^{m) — d^{m,k). We have thus determined the 



counterterms when m is half- integer. 

Notice that when m is half-integer, a finite counterterm is also possible: we may intro- 
duce a nonzero am-i/2) producing a finite counterterm proportional to c^i^(— A:^)™"^/^ c^. 



Something similar happens for scalars [|61| . The scalar equation of motion has Bessel- function 
solutions, and for certain values of the scalar mass m the linearly-independent solutions are 
Bessel functions of the first and second kind, which means the asymptotic expansion includes 
logarithmic terms. In those finite counterterm also appears, and is proportional to 

the matter conformal anomaly of the dual CFT [^] (and references therein). Here we are 
simply seeing the fermionic version of the scalar story. 

To summarize: for the action in eq. ( p. 2D , clearly the bulk term vanishes when evaluated 
on a solution. The only nonzero contribution to the on-shell action comes from the boundary 
terms, hence S = S^dy on-shell. Generically, Sbdy diverges and we must add counterterms. 
When m is not half- integer, the resulting on-shell action is^* 

S = J d'^x cfc^ + O (e^) (A.38) 
while when m is half-integer but m ^ 1/2 (here we take d^im, k) — )• d^), 

S = I d'^xc^c^ (^-i+^^-i+E^^^^;.-!-.) {-kT"'+0{e'lne), (A.39) 
and for m = 1/2, 

S = j d'^x C+C2 + O (e^ In e) . (A.40) 

These on-shell actions remain finite as e — )• 0, and hence, upon functional differentiation, will 
produce renormalized field theory correlators, as we next discuss. 

A. 3 Computing Renormalized Correlators 

The field \I' in the bulk is dual to some fermionic operator O. As reviewed in section ^, 
the renormalized on-shell bulk action, Sj-en = li™e-s>o 5", acts as the generating functional for 



''^Recall once again that we restricted to positive m. At the end of the last subsection we explained how to 
recover results for negative m from the results for positive m. 
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correlators involving O. In other words, to compute renormalized correlators of O, we take 
functional derivatives of S with respect to some source. We identify the source for O as the 



coefficient of the dominant term in ^''s near-boundary expansion. In eqs. ( A. 11 ) and ( A. 15 ), 
the dominant term is the u2~"^ term, hence we identify C2 (fc) as the source for O. More 
formally, we equate (see eq. ( 2.20 )) 



j (fx (c^O + Oc^) 



(A.41) 



where the left-hand-side in the exponential of minus the on-shell bulk action in eq. (|2.2| ), and 
the right-hand-side in the generating functional of the dual field theory, with c^(fc) acting as 
the source for the operator O. Upon taking minus the logarithm of both sides, we find that 
the on-shell bulk action is the generator of connected correlators. 

For any value of the bulk fermion's mass m, a quantization exists for which O has 
dimension A = ^ + \m\ |Tl|, When \m\ G [0, 1/2), a second quantization exists for which 

T 



A = i — \m\. These two quantizations correspond to two different field theories, one of which 



(the one with A = ^ — \m\) is "unstable" against a relevant deformation by a double-trace 



operator of O, and flows to the other, "stable," theory, as in the scalar case |104| When 
\m\ > d/2, the dual operator is irrelevant: A > d. 

When d is even, O will be an operator of definite chirality, since in that case 7" is the 
chirality operator of the field theory and c^(A;) has definite chirality. Notice that taking 
m —7- —m takes c^(A;) — t- cf{k) when m is not half-integer and c^(A;) — s- (— l)l™'l'*"2 c^ik) 
when m is half-integer, and hence switches the chirality of O. 

We reviewed in section ^ how to compute the Euclidean Green's function from the renor- 
malized on-shell action. In particular, from eqs. ( p. 251) and (|2.26| ), we have 

(00) =-^^^ = -^ = G{k)j\ (A.42) 



where G{k) is defined by cf = — G(/c)7*C2 . As we mentioned below eq. ( 2.26| ), in general. 



we extract G{k) 7* from from a solution by imposing some regularity condition in the bulk 
of the spacetime. As an example, let us consider the simplest case: pure AdS^+i with m 
positive and not half-integer. The solution in eq. ( [A.5| ) is then the solution for all u, not just 
the asymptotic solution as u — )• 0. We need to write cf in terms of C2 ■ To do so, we impose 
a regularity condition in the bulk of AdSd+i- Deep in the interior of AdS^+i, where w — )• 00, 
the solution in eq. ( [A.5| ) diverges unless = (— l)™"^^/^ (7^. Translating that condition into 
a condition on and is trivial, once we recall the definition of and in eq. ( A. 91 ). 
We find 

c+ = -2-2™ (2m - 1) ^ ~ ""j A:^— 1 c+. (A.43) 
r (2 + m) 



''^For the Dp-branes we study in section 4.1, m is always integer or half-integer, so the second quantization 



is only possible when m = 0. Indeed, in the bulk of the paper we use m = 0. In that case, however, the two 
quantizations are equivalent (see appendix A of ref. and ref. |6c|]). 
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We then use eq. ( A.13|) to write in terms of i^C2 . We finally obtain for the two-point 
function, 



[OO) 



-2m 



r(^ 



k 



2m- 1 



(A.44) 



r(i + m) 

which agrees with the result in refs. |62, |6^, |105(| and, up to normalization, is the correct 



momentum-space form for the two-point function of a quasi-primary operator of dimension 
A = d/2 + m. Similar arguments work when m is half-integer, where again we find the correct 
momentum-space form for a quasi-primary operator of dimension A = d/2 + ra. In that case, 
we can also adjust the normalization to any value we like using the finite counterterm. 
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